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Preface 


It has been noted, in the course of study of geodesy, that 
there is a sparsity of material devoted to the mechanics of determining 
the size and shape of the reference ellipsoid from actual observations. 
Bare outlines of some methods are given in the few current textbooks 
available in English, and although some of the older texts, notably 
Clarke (5) and Crandall (6), do give examples, the methods used are 
rather ancient and OU mg ana in notation systems quite unfamiliar 
to the present-day student.  Hayford's original study of 1909 is the 
best available (10). Since one of the basic aims of geodesy is to 
determine the figure of the earth, it seemed like & good idea to 
pursue the matter, if even yet somewhat superficially, to the point 
where there would be few questions left concerning the manner of going 
about the task. 

The writer wishes to acknowledge the help and advice provided 
him in preparation of this thesis, particularly by his adviser, 

Dr. W.A. Heiskanen, Director of the Institute of Geodesy, Photo- 
grammetry, and Cartography, the Ohio State University. Acknowledgement 
is also due Dr. U.A. Uotila and Mr. Ivan Mueller, both of the Ohio 
State University, for certain helpful suggestions; the U.S. Army Map 
Service and Mr. W.M. Kaula, for data on points in the Finland - 

South Africa arc; the U.S. Navy, without the instigation of which, 

the course of study referred to in the opening sentence would not have 
been possible. 
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Finally, the writer wishes to express his appreciation to 
Frauke and Peter Wilson, the former for her fine job of typing the 
bulk of this thesis on short notice, and the latter for certain 


editorial assistance to the typist. 
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Chapter 1. Introduction 


1.1 Preleminary Remarks. 


That the most convenient practical reference surface for purposes 
of geodetic computation is the ellipsoid of revolution has long been 
recognized. This is emphasized by the fact that the major datums in 
use in the world today are all based on ellipsoids of revolution. The 
fact that this shape is a relatively simple regular mathematical sur- 
face, and that it closely approximates the mean geoidal form, make its 
choice quite understandable. The Russians, who have been proponents of 
a triaxial ellipsoid, do not even consider its use for practical 


geodetic computations. Izotov (18) states:- 


"The triaxial ellipsoid could not be accepted for triangulation 
processing, as this would have made the working out of geodetic 
measurements too complicated”. 


To dismiss the subject of triaxiality for purposes of this paper, the 
following is cited from "The Earth and its Gravity Field", by Heiskanen 


and Vening Meinesz:- 


"It appears that we had better forget the triaxial ellipsoid, 
particularly because since the earth is in close isostatic 
equilibrium, it can hardly have the shape of & triaxial 
ellipsoid." (13) p. 79-80. 


The use of geodetic data to determine relatively precisely what 
size and how much flattened the reference ellipsoid should be, is the 
subject of this discussion. In this connection, one of the principal 


purposes of geodesy has been variously stated: - 
I 





"Its theoretical function is to determine the size and shape 
of the earth ..." (13); "Geodesy is the science which treats 
of investigations of the form and dimensions of the earth's 
surface." (17); "The accurate determination of the figure of 
the earth." (6); "In combination with observations for latitude, 
longitude, and gravity to assist in determining tbe size and 
shape of the earth, ..." (2); "It may be said that the ultimate 
aim of scientific geodesy is to determine the size and shape of 
Che geoid, ..." (2). 


Of course, data collected in geodetic surveys along with the calcula- 
tions connected therewith, and computations and results of computations 
of the figure of the earth are mutually interdependent. 

Although the problem of determining the best size and shape 
for the reference ellipsoid has decreased in importance relative to 
other modern forms of geodetic endeavor, there are still significant 
efforts being made in this field, notably by the U.S. Army Map Service. 
This is not to say that the exact dimensions of the earth are not im- 
portant, but the position has now been reached where further precision 
given to a reference ellipsoid must be subordinated to the problem of 
determining the deviations of the actual form of the geoid with respect 
to the reference surface selected. 

This change of emphasis is described somewhat succinctly by 


Tengstrom: - 


"The principal problem of geodesy is to determine the geocentric 
coordinates of all points on the earth's surface, a task which 
was earlier referred to as a determination of the 'size and 
shape! of the earth." (30); 


or by L.G. Simmons:- 


"The ultimate goal of the geodesist is, I suppose, to determine 
the parameters of an ellipsoid of revolution which best fits the 
figure of the earth as a whole. But this is not all. He also 





is concerned with the details of lack of fit of this ellipsoid 
to the actual Earth's shape, the geoid." (28). 


Most important would be the selection of one reference surface 
to be used in problems of world-wide extent, such as for requirements 
of missile warfare, the success of which might very well be determined 
by accurate knowledge of locations of launch sites and targets.  Natu- 
rally, the best reference surface obtainable is to be preferred. The 
establishment of a "World Geodetic System" and conversion of existing 
systems to it, has for several years been one of several significant 
projects promoted and encouraged by Dr. Heiskanen (16) (13) p. 299-310. 
Both the U.S. Air Force and the U.S. Army Map Service are developing 
world geodetic systems for military purposes. These are eventually to 
be combined into a single unified Department of Defense World Datum, 
the contributions to the establishment of which will be made by all 
U.S. military services. The reference ellipsoid to be used for this 
will be based on all available information from all parts of the world. 
(3) p- 77. 

Having exposed the near-obsolescent character of the problem 
of determining the dimensions of the reference ellipsoid, we may now 


proceed to describe the several variations of its solution. 


1.2 General Discussion of Methods. 


The standard approach to the problem under consideration has 
been, from earliest times right to the present day, that contained in 
the arc method. Variations of this approach have become refined and 


quite sophisticated, taking into account the earth's isostatic 





equilibrium, and in some cases being combined with the gravimetric 
method. The earliest recorded application of the arc method is told 
in the now trite story of Eratosthenes' activities in this field. 

The arc method assumes that short meridional arcs on the 
earth's surface are arcs of circles; when two meridian arcs are 
measured, one at relatively high latitude, and one near the equator, 
both the equatorial radius and the flattening can be computed. The 
more such arcs used, combined in a least squares solution, the better 
accuracy is obtained. 

A more convenient adaptation of the arc method has been used 
since early in the present century. This uses a combination of 
geodetic and astronomic data in the form of astro-geodetic deflections 
of the vertical. Hayford used this method in developing what is now 
known as the International Ellipsoid (10). He called the procedure he 
used the "area method”, because the data was taken from many points in 
a large network of triangulation. He made several solutions, the 
deflections of the vertical in each being reduced to the geoid (co- 
geoid) according to the theory of isostasy using different depths of 
compensation. He chose the solution in which the sum of the squares 
of the residual deflections after adjustment was the least. Thus a 
"by-product" of his work was a solution of the mean depth of compen- 
sation for the earth's crust. 

The application of the gravimetric method to the problem under 
consideration is really a combination of the arc method, modified by 


the introduction of gravimetrically determined deflections of the 





vertical, to obtain the earth's equatorial radius, and the deter- 
mination of the earth's flattening from gravity anomalies by use of 
Clairaut's formula. 

The latest, and of course the most modern method used to 
determine the earth's shape, is the artificial earth satellite method. 
By measuring the perturbations in the satellite's orbit around the 
earth, caused by deviations from the spherical in the earth's shape, 
the flattening is obtained. „This method is so new in actual 
application, that an evaluation of results thus far obtained can only 
be the subject of conjecture. 

One other method, similar to the deflection of the vertical 
application of the arc method, involves the use of the geoid un- 
dulations, obtained both gravimetrically and astro-geodetically. The 
solution is obtained in reducing the sum of the squares of the 


differences between the two at many points to a minimun. 





Chapter 2. Dimensions from Two Arcs 


2.1 Meridional Arcs. 


In this analysis, as in all subsequent ones in this paper, the 
reference surface is assumed to be an ellipsoid of revolution. The 
mathematical properties of an ellipse form the basis for expressions 
referring to the meridian of such a surface. 

Tne measurement of two arcs in meridian, one near the pole and 
the other near the equator, eds the data for this, the simplest 
means whereby the earth dimensions may be obtained. The arcs, which 
may be on different meridians, are measured, and the difference in 


latitudes between extremities for each is computed, the latitude of 


one end point for each arc having been determined astronomically. 





Meridian 


Fig. 1 
In figa L, 94 and Po are the mid-latitudes, and $1 and S, the 
measured lengths of the two arcs. Assuming S, and S, to be arcs of 


circles, we can write:- 





ell 
2-1 R- = —— ese 1" 
S 
2 
R. = — cse l” 
2 An 


AP: and AY, being expressed in seconds. 
The expressions for the radii of curvature of the elliptical 


meridian at the respective latitudes are: 


SERES 


ae) Ryn — rO 
(1-e“sin P) 
R Mies) 
2 3/2 


(1-e°sin o) 


Equating the two expressions (2-1) ana (2-2), we obtain: 





S 2 
(2-3) A esc 1" = — 
1 (1-e “sin P) 
le esc 1" = alle) _ 
bà 3 
Mo (1-e"sin“9,) “ 
Dividing Ry by R, in the form (2-3) and solving for aer 
S. Mp 1-e“sin“q 
ID 2 a 2 
E) = 7 ` 4 
2 1 l-e sin P 
e 


ae " 159 
sinp -q sin p, 


Then, the value of a can be found from (2-3): 
2.2 %2 2.2 1% 
s, (1-e sin 0. ) s (1-e sin Po) 


(8-5 ^ & » ——M——— -= > 
Ap, (1-e )sin 1" Ap (1-e lisis! 





Now all that remains is to determine the flattening, which is 





simply: 
a-b b 
(2-6) E clo x ERE 


2.20 Obliqme Arcs. 

In the early efforts of measuring the earth's figure relatively 
short arcs were measured, almost always in meridian. In later deter- 
minations, existing triangulation data was put into use, and consequent- 
ly successive points used were not always on the same meridian. Such 
geodetically computed arcs, chosen so as to be inclined to the meridian 
at as small an angle as possible, could be reduced to the meridian by 
appropriate formulae or from tables based on the used ellipsoid. Thus 
armed with the same data as in section 2.1 a similar problem can be 
solved. 

When the length of the oblique arc is measured, along with its 


azimuth and mid-latitude, its meridian distance is given by: 


(2-7) MUELsin 1! = - s cos A - am g^ sin“A tan Q 


+— 83 sinóa cos A (1+3tan o) 

mo 

ÓN 
where N is the radius of curvature of the prime vertical and A is the 
azimuth, reckoned positive clockwise from south. 

Using (2-7) a chain of triang.es can be projected to the meridian 

for the whole length of the chain. Given geodetic latitudes, However, 
use can be made of specially prepared table, such as (34), (35), or 


(36), which have meridian distance tabulated against latitude for the 


ellipsoid involved. 





J 


The above method of determining the earth's dimensions (Sections 
2.1 and 2.2), though it represents the basis for other methods, is 
mainly of historical and mathematical interest, since no longer are 
serious scientific computations carried out in such simple terms. 
References for the presentation given in these two sections are (6), 


(17), and (23). 





Chapter 3. Dimensions from Several Arcs. 


3.1 General. 


When geodesists of the 18th Centurv first began to compare the 
results of their arc measuring efforts, they found such large discrep- 
ancies that they began to speculate about tne accuracy of their meas- 
urements and &bout the correctness of the assumption that the earth's 
true shape is an ellipsoid of revolution. Of course the accuracy of 
surveys of those early days was not what it is today. Nor couid one. 
hope to obtain the same results from arcs, taken two by two in dif- 
ferent parts of the world having varying kinds of topography, meas- 
ured by the crude methods of the day. Not only was the measurement of 
geodetic triangles so rude that spherical excess remained undetected, 
but the process of adjustment by the method of least squares was un- 
known (23). 

Obviously, not wanting to discard such a convenient reference 
surface as an ellipsoid, which after all departs from the geoid sur- 
face by probably less than 100 m. at the maximum, and from the 
equilibrium figure by even less (of the order of 2m. at = 0,45, 
and 90, (13) p. 59), the best thing to do would be to combine the 
results of several arc-measurements in such a way as to determine the 


mean values for the earth as a whole. 
3.2 Laplace's Method. 


Toward the end of the 18th Century, Laplace published the re- 


sults of his attempt at combining the data from nine different sources, 


10 





Ja. 


extending from Lapland to Cape of Good Hope in latitude. His method 
involved the data corresponding to an arc of one degree of latitude for 
each of nine different meridians. 

Again considering these one degree arcs as arcs of circles and 
using R, the radius of curvature of the midpoint of the arc, and 9, 


the mid-latitude, we have:- 
2 
2 71 R T & (1-e 
bay as BPR o zeae’) 
910.1. sọ” 
Developing this in series, we obtain:- 


e 
(3-2) a- Tele) (+20 sin? gi Be! sin' o+...) 


which can be written in the following form:- 


(3-3) dei ai) sins Dar E Xs D + eee 


Tra (o5) 


in which M = A 
a 
N = 5 e M 


P= 


l 
(b 


M, o... etc. 


Powers of e higher than the square are neglected. 
In order to combine the data for each degree -measurement in 
the manner yielding the least discrepancies, the observation equations 


for each arc were of the form$i- 


(3-4) d = M-N sin D 


AA 
d, -M =- N sin 9= V2 etc. 


in which Vi» Va» etc., are the residual errors. Since the method of 


least squares was unknown at the time, Laplace combined the error 





He 


equations in such a way that the absolute values of the residual errors 
would be a minimum, and the algebraic sum equal to zero. Then solving 
simultaneous equations for M and N, the values of ef and a are found. 
Of course, this problem is easily adapted to least squares solution. 


The results of Laplace's computation are:- 


M = 11 061.2 meters 
N = 1 194.9 meters 
ES = 0.007262 

a = 6 383 600 meters 
f = 1/278 


The error, or difference between a measured degree and one 
computed from the result of the crude adjustment of the nine arcs, is 
of the order of 268 meters. From this, Laplace concluded that the 


earth deviates materially from an elliptical figure (23). 


3.3 By Correction of Latitude. 


A refinement of Laplace's method, in use by the end of the 19th 
Century, is here described (6), (17). This also involves the formation 
of observation equations between observed latitudes and either directly 
measured meridian arcs or oblique arcs projected to the meridian. 

For short arcs the relation s = Ap R sin 1" can be used, but 
for longer arcs a correction is required. Crandall (6) states that the 
simple relation given is good for a Ap of "several degrees" because of 
the large probable error of latitude determination, some O?0h (more 


likely this will be on the order of 01 at best), in addition to the 
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error caused by local deflections of the vertical. The correction, for 
use with longer arcs is:- 
3 


(3-5) as = ga e^ (apt sin 1")" cos 2 p 


Given astronomic latitudes for several stations in different 
arcs, as well as geodetically measured segment lengths s, we choose 
separate initial points for each arc. From the first two latitudes in 


any one arc we have:- 


S d s 
(3-0) 9, - 9 * goi 1" ^ Fin T^ 
where 
u e^ sin* 3/2 
(3-7) IB = 2 
a (1 - e^) 


In this method we make use of the used reference ellipsoid, for 
which we know the values a and es, and compute corrections da and 


Se“, to that ellipsoid resulting from the data given. Thus a = a + Sa, 


Expanding equation (3-7) into a series anà substituting it in 


equation (3-5):- 


2 
(3-0) Da T Re ee + (1 - 5 sing) | A 
O 


in which Ro is the value of R for a and eS 


Considering the measured meridional arcs accurate in comparison 
with the observed latitudes, we give corrections E V5» to latitudes 
Q4 5 Pos etc. Since in the derivation of (3-8), all terms containing ` 

2 2 


were neglected in the coefficients of Da and Se”, i.e., E7 was con- 


sidered equal to zero, we can say R = a and s/(a sin 1") = Po = P 
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Thus equation (3-8) becomes: - 


D S S Sa 3ó 2} 
(potv,) - (p,+v7) = R sin T" E & Sin 1" E a +(1- ¿sin q)Be 


d s 
- R. sin 1" ; and then 


o E WA - 89 a cae 
2 1 a sin 1 2 2 l 
S - ds 
ii R- sin 1" ^ (05 - 94). 


* 


This becomes the basic observation equation to be written for 
each segment of each &rc. The unknowns are ME 5 a, and 5 E. 

In the least squares solution, in which ENS - minimum, the 
magnitudes of the various coefficient elements make it convenient to 


use the following modifications. 


Placing X -65 aj 1000 y = 1000 8 e 
m, = - 1000 (9, - 9,)/a, 
9, - O 
JM E 3 2 
n, 1600 (1 5 sin p) 
ie ei de a lo - ) 
2" R sin Da 


We can then write a simple error equation for each segment, 


substituting the above symbols in (3-9):- 
vv tm xX+nm,yt+t l5 
Since the initial latitude must be corrected, an additional 


1 = Pte In this 


method, there will be as many unknowns as there are arcs plus two, i-e., 


observation equation is written for it, simply v 


the v for the initial latitude of each arc, plus x and y. In the case 
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that a value for the initial v is assumed to be some arbitrary value, 
e.g. zero, there will be only two unknowns. 


For each arc, observation equations will be as follow:- 


(3-10) Vi zum 
Va + m,x + ny + l, = Ji 
+ mx + ny + 1, = 
NS E i 
X XA x x x 
v, + mx + ny t ` „`? 


The result of a computation of this type presented by Crandall 
(6), and using data (circa 1890) given in the “Handbuch der Vermessungs- 
kunde" (37) for six arcs (French, English, Hanover, Prussian, Russian, 


and Swedish), varying in length from 095 to 127%, is as follows: - 


a = 6 377 800 meters 
ae 0. 0069091 
f - 1/289 


The standard error of the mean amounts to 271 for & latitude 
determination, which is quite a bit greater than the accuracy obtain- 
able for first order work, (about O71). Thus, Crandall concludes (6), 
an ellipsoid of revolution will not fit the data without large local 


deviations of the vertical. 


3.4 By Correction of Radius. 


A method given by Bomford (2) p. 364, which is a variation of 
that described in Section 3.3, works with radii of curvature at points 


in an arc rather than with latitudes. 
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For given values of R in widely separated latitudes, the values 


of a and e are obtained from the expansion of equation (2-2):- 


2 
3-11 = a (l - e”) 
EL (1 - e sin“ ọ)3/2 


a (1 - em) (1 + 2 ef + de ii panas) 


Using a =e fe £^ and the approximation 9 ae eot, and 


droppiag the terms of higher order than gen the following form is 


obtained: - 


* 


(3-12) oR-628-28 6f (1 - 2 sin^q) 


in which 6 R = R - Ro results in changes 5 a and 5 f in a and f f 
the ellipsoid provisionally adopted. From equation (2-3) we may use 
R, = $ csc l"/Ap' for ares shorter than about 100 miles (2). The 
value of Ro would be obtained for the used ellipsoid. 

In a more practical application, more than two values of R 
would be available, and observation equations for least squares 


solution would be as follow:- 


3 2 s,csc 1" 
x = ea y (1 = 5 sin qi ) - man T (RAN 


ll 


(3-13) v 


s cse 1" 


e 
Po ) d a. + (Roo 


< 
il 


3 
eS (1 - 5 sin 


x ¥ He x x 
s ese 1" 


Ie 
X - 28, y (1 - 5 sin 9, ) - T * us 


ll 


V 
n 


where x = O6 a and y = ô f are the unknowns, s. is the measured length 


l 
of the arc P4P5; and Ap" = Po = P in seconds, etc. Again X v^ is 


made & minimum. Weights may be assigned &s deemed necessary. 





Chapter 4. Dimensions from Arcs, Using Deflections 


Greeule Ver unea| 


4.1 Basic Considerations - Astro-Geodetic Deflections. 

(Note: A more complete analysis of "basic considerations," in- 
volving absolute deflections of the vertical, is given in section 4.3.) 

As mentioned earlier, a very convenient adaptation of the arc 
method makes use of astro-geodetic deflections of the vertical, which 
are generally available in profusion wherever triangulation control nets 
have been made. This has been done in many parts of the world in 
connection with the various Onell surveys being Carried out. 

This method can be applied to individual arcs of meridians in 
one extreme, and to the whole network of triangulation points in the 
area covered by a national control survey in the other. Because of the 


"area 


latter application of the method, Hayford chose to call it the 
method" when the results of the U.S.C. & G.S. computations of the figure 
of the earth based on data collected in the large area covered by the 
United States were published (10). 

Here again, the reference ellipsoid used in computation of the 
geodetic data is used as a point of departure, or more accurately, an 
"area of departure", from which we must apply measurements and compu- 
tations to get as close to the best ellipsoid as we can. This "best" 
ellipsoid is the one most nearly approximating the earth spheroid or 
mean geoid surface, in size, shape, and orientation. 


Figure 2 shows in simple terms, the relationship which can exist 


between a reference ellipsoid and the geoid. 


1 
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Fig. 2 - Astro-Geodetic Deflections of the Vertical 


In general, thereware six diiterent items Chat Way be corrected 
by the procedure involved, and the particular ones that are corrected 
depend upon what is desired and upon what and how much data is supplied. 
The existing geodetic data, i.e., latitudes, longitudes, deflections, 
etc. ¡depend upon the initial point (datum) of the survey, and upon the 
assumptions made about that point. The six items are: the equatorial 
radius a; the polar flattening f; and the following quantities assumed 
forthe initial point: its l8 1ï 9 its longitude N3 the geoid 
height Nos and the azimuth A of its meridian. In many cases, all that 
are desired would be the corrections Sa and óf. Usually the cor- 
rections to Po and Ne at the initial point will be in the form of 


corrections to the meridian and prime vertical Components to the 





Bm 
deflections of the vertical. 
The following notation will be used in the discussion to follow: 


Qo - Astronomie latitude and longitude 


Q, A Geodetic atitude and longitude 
(o! -q) - Meridian component, astro-geod. defl. 


(A' -A)cosp - Prime vertical comp., astro-geod. defl. 


9 - Residual deflection in meridian 
») - Residual deflection in prime vertical 
đến - orya gaien corrections to deflections 
E due to displacement of ellipsoid 
2 
đo - Change of ellipsoid corrections to deflec- 
7 tions. These contain the values ôa 


and of which result when an ellipsoid 
is changed. These corrections actually 
apply to the geodetic parts of the astro- 


geodetic deflections. 


When it is desired to find a better orientation of the ellipsoid 
as well as to find its size and shape, the observation equations which 


may be written for each point take the following general form: 


(4-1) € 


li 


(e -9 *8$5 * a£, 


i] 


/ 


3 
(A - A) cos P+ 47) +47, 


If it is desired only to find a and f, using data from a meridian 


arc only, the series of observation equations will be as follows: 
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(4-2) E, = (o - o)+ (4 8,), 


vn 
PO 
N 
na 
-9 
PO 
t 


Pat (4 57), 


2e 96 ee e 


uve 
n 
A 
O 
Ha) 
B 


A 


Then the least squares solution will involve making zee = 
minimum, this giving us the parameters a and f of an ellipsoid which 
will best fit the condition which, in effect, forces the residual de- 
flections to a minimum, as ose to zero as possible. Weights can hs 
assigned to the observation equations as deemsd necessary. 

In general form, the condition can be stated: 

(4-3) AC. + y ^) - minimum, 
which is good for any problem. Either or both terms, E^ and »*, can 
be used according to the problem to be solved. 

In this case, in which astro-geodetic deflections are used, the 


better results are obtained when these deflections have been reduced 


isostatically. 


4.2 Formulas for Displacement and Change of Ellipsoid Corrections. 
There are four well-known methods of computation of the change 

of ellipsoid and displacement corrections to the astro-geodetic 

deflections. These are based on formulas representing different 

derivation approaches to the problem. They are the formulas of Hayford 

(10), the U. S. Army Map Service (+), Bomford (2), and Vening Meinesz 

(33). The derivations of these formulas are given in the references 


cited (except (4)), and being beyond the scope of this paper, will not 
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be covered: here. 
4.21. Hayford's Formulas. The form of otservation equation given 


in (10) is as follows: 


(4-4) =k (p) +1 (A) +m (A) * nj (485) * o, (10,000 e^) + (4'-9) 
4 =k, (9) +1, (1) + m, (A) + 0,48) + 0, (10,000 e”) 


ua! Sb coB 9 


in which (q), (A), (A) &re unknown corrections to the initial q,A, and 


a 


100) and (10,000 eS) are unknown 


A at the origin of the survey, and ( 
corrections to the values of a and a for the reference ellipsoid. It 
can be seen that these are similar in form to equations (4-1). 


The change of ellipsoid corrections are: 


a 2 
(4-5) dga = n (o + 0 (10,000 e”) 





a 2 
dy. = Ny (150! + 0, (10,000 e”) 
i : _ 100 aN 
in which Ra = n ae 6 COS Az 
2 e 
S e 6 2 2 
8 = — |1 + ——s— COS Q_ COS 
N | 6(1-e°) O Ân 
sin Q9 
N O 
A TA 8 cos 
+ - 20,000 sin 1" (1-e sinó) = | 
R ph C 
Y] AAA (o-q_)[1+2 sin p_ + 3 cos( p +9) ] 
40,000 sin 1" (l-e`) 
100 : 
d we. ccm = 
sino, 
OLE —l [6 sin 
S 20,000 sin 1"(1-e sin) E 
(1-e2sinp )1/2 coso, 3 5 
+ 907 cos Ap sin A, 


60,000 a sin 1" (igen 
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The corrections for displacement (orientation) are:- 


(4-6) — aty-k, (9) «14 (A) +m, (A) 
an, = Ea (p) + L, (X) * m, (A) 
Pa 
_ sin A) 
in which ki = =- l + E (1 + Q) `. TR 
2 sin 5 (A, > An) 
l = Zero 
sin [1 + cos (A - A )] 
m = en 1+0 A eee 
F R 2 sinf a - A.) 
5 Vg - ấn 
k, = - sin 9 sin SA - A) 
l, = = cos (0 


cos @ cos A, sin (X - A) 
2 sin Ar 
ge 


Ecos? lay ap) 


3 
i 


D 
N 


The purpose of the coefficients k, 1, m, n, and o can be stated 
as in this example for one of them: ki is & numerical coefficient, 
such that if the latitude of the initial point were corrected by the 
amount (qp), the change produced in q' - q would be ! (op). 


It can be seen from inspection of the coefficients that for 


points on the same meridian as the initial point, the values are 


kl = -1, 1, = O, and m = O, so that the observation equation would be:- 
E dl = i 
(4-7) E m" (p) + aT (356! + eF (10, 000 e ) T (9 = 0). 


4,22 Army Map Service Formulas. In an article in Transactions of the 


American Geophysical Union, No. 37, 1956, Chovitz and Fischer 
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reported on the U, S. Army Map Service's determination of the figure of 


the earth from arcs (^). The procedure is described in some detail, 
the essential part regarding correction formulas being reproduced in 
somewhat different form here. The method is designed to be correct for 
meridional arc segments of r but not more than cu in length, and for 
arcs of parallel of comparable length. 


For meridian ares the formulas are: 


(4-7) du = En - 9) + de [1(9,-0,) i (9$ - 94.3) 


* 


[11 cos 20. - III cos AQ + IV cos sonl} ese 1" 


mid-latitude of segment = 


bar 


in which q, (9 4 + P) 


(4-8) at 


- ©.) = initial orientation of the arce segment. 


il 
vr 


D 


In this method the computation proceeds segment by segment, with 
a new initial point for each segment of the arc when the length limits 
given above are reached. Thus the procedure is similar to the one 
described in section 3.3. This presupposes that all the points in the 
arc furnish data that is adjusted to the same datum. In contrast, the 
Hayford formulas (and those to be given subsequently) are designed to 
cover a large, but limited area, using only one initial point in the 
computation. The unknowns in formulas (4-7) and (4-8) are ES» da, and 


Dess The Roman numeral coefficients are: 


aj l 4 

I m T. e^ - zí e hee = ica 
MP s prpe “ae er 
3 d 9 4 
SEDE = He gee 
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For arcs of parallel the formulas are: 





$a 2 sin“ p cos Q 
(4-9) đe = | 5 AA) cos 9 + de (RR) er E e 
Gace sin^g + e sin'g)| eso 1" 
where © = latitude of the parallel. 


ESA uu 9 
(4-10) ad 13 "fs (N ho) cosp - initial orientation of the arc segment 


The unknowns are N >» 98, and ss Solutions for meridians and 


parallels can be done separately or can be combined for a cross-solution 
(4). 


4.23. Bomford's Formulas. The formulas for change of ellipsoid 
given in Bomford (2) p. 127-130, actually contain corrections for both 
displacement and change of ellipsoid. They are developed here in 
separated form. These formulas are used in Bomford's version of the 
area method, which supposes a “homogeneously computed national survey" 
covering a 20° by 20° area containing 50 astro-geodetic stations. 
Bomford states that it will generally be best to ignore values of X% 
deduced from azimuth observations (these were omitted in the description 
of Hayford's method above), and to use only those based on longitude 
(2) p. 365. 

The change of ellipsoid formulas as given by Bomford (2) p. 129, 


365, are as follows: 


(4-11) de =- | (1/a) six v cos w (U sin u + V cos u) 


+ sin u sin 0 óz. - (1/a)cos u(V sin u - U cos u_) 


+ sin 2u se| ese 1" 





eo 


- = = ] i = tí 
(4-12) dy | (1/2) sin o (U sin u, + V cos u) cos w 59, | ese 1 
where U = -a SE - aóf sin 2 ul 

ee & 

O O 

u - reduced latitude; tan u = (1 - f) tan 9 

® = longitude measured east from origin 

N = height of origin in corrected ellipsoid 


above that in reference ellipsoid 


ôg and ön are in radians. 
O O 


Substituting the values of U and V in (4-11) and (4-12), and 


separating the expressions into two parts: d&g = d&, + d&, na 


d» = dX. + AND we obtain: 


| 


(sin u cos wm cos u, - Cos u sin u) daja 


(4-13) de, 


+ [sin u cos w sin u cos u + cos u sin u 


(1 + cos^u ) - sin 2u | d ese 1" 


(4-14) d£, 


(sin u cos o sin u, 4 cos u cos u ) óE " 
O O O 
- sin u sin woy"- (sin u cos w cos u 
O O 
- cos u Sin u ) (N /a) ese 1" 


: e "i 
| sin w cos u Sa/a) + sin w sin u, cos u, sr] ese 1" 


(4-15) d No 


(4-16) d%, = sin o sin u SE, + cos w on, 


- sin w cos u (N /a) esc 1" 


4.24, Vening Meinesz' Formulas. The formulas here referred 


to are those published in 1950 (33), and which were developed for the 
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purposes of (a) describing the relationship between the variation of 
the vertical deflection components at the origin of a survey and the 
changes resulting in the deflections at any arbitrary point in the 
survey, the variation having been caused by a parallel displacement 
of the reference ellipsoid with respect to the geoid; and (b) 
describing the changes in the deflections of the vertical brought 
about by a change of reference ellipsoid. 


The complete formulas referred to in (a) above, are as 


follow: - 
Ro 21 
(4-17) ee = p iun) 33 sin 9 sin c, sin 5(A-A, 3) SE 
O 
+ =~ sin 9 sin (A-A,) UM 
T ; : mod 
* E [sin (9-9. ) - 2 sin 9 cos q, sin 5 (A-A JONS 
To 
(4-18) an, =- sin 9, sin ( A-A) * — cos AA) en 


le n] 0 
O 


cos 9, sin ( A-Aa) SN, 


where R, R = radius of curvature in meridian at initial point P_ 


and at point P, respectively 


5S 
= 
II 


E radius of curvature in prime vertical at initial 

2 

point P, and at point P, respectively (r is used 
here, rather than the usual N, which is reserved 


for geoid distance). 


The above form is that given by Heiskanen in (13) p. 308, and 
is better suited to the use of tables than the form in the original 


paper (33), which uses all-e? J/W rather than R, and a/W rather than 





eT 


` 


r, where W = esto)”. 

Vening Meinesz points out that his formulas are similar to 
those developed by Helmert (37), but include SN, terms as essential 
differences. 

If the line POP is long, the complete formulas above must be 
used. However, for nets of small east-west extent, terms 
proportional to the squares of (A-A,) can be neglected; and for 
nets of small extent in latitude, the ratios R/R» r/R, R /7, and 
pt are very nearly 1 (within 0.01), and the formilas (4-17) and 


(4-18) become: - 


(4-19) de, = cos (0g BE * sin @ sin (AA IN 
O 
+ sin (9-9. ) "m 
(4-20) dn, - - sin 9, sin AR) BE + 8n. 


ON, 
- cos 9, sin AM) A 


According to (13), the formulas (4-19) and (4-20) apply when 
pS 300 km. in east-west direction and < 1500 km. in north-south 
direction. However, Vening Meinesz (33) states that if it is 
desired not to neglect quantities greater than 1/300 of the 
deflection of the vertical, these formulas can be used for nets 
extending up to 4° ho' in east-west direction and to more than 20° in 
latitude, from the origin ps Furthermore, the complete formulas 
will seldom require that the ratios between R, Ro? r, and ro be any- 


thing but 1 even for large nets, to attain the required accuracy. 
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The complete formulas referred to in (b) above, for change of 
ellipsoid, are here transcribed: - 


(4-21) de, = (Esto (9-9) - 2 cos p_ sin y sin” > (A )] op 


l ar 
- 4 cos y cos z (ptp) sin 5 (p-p) Sf 


+ 


. 2 . 
[(2-3 sin p - cos (pq+p_)) sin (9-9,) 
€ Pr o 2 2 o 2 1 2 ` 

- ? cos q, sin @ (2+sin @_ - 3 sin o)sin 5 (X A) 

o o 2 ji l f 
- sin 29. sin” 5 (9-9. ) - 2 cos @ cos 5 (9*9) 

o 1 1 a o 2 l 4 

sín 5 (9-93) f5B + [ 5 sin 4 9c (sin > 9-9.) 

o o e 2 i o 3 
+ sin 9, sin sin = (NA, )) - sin 9 sin > (99) 

B li mm 
sin 5 (9-9. ) - sin 9 sin > (op) sın 5 9-9.) 
- 4 cos y cos L loro ) sin 2 (9-9 ) 
2 O 2 O 


2 1 xc " 
+ 4 sin 2 y cos > (pP) sin 5 (9-9, )118£ ) escl 


(4-22) an. = {- cos q sin (A-A ) 8Ð 
+ cos 9, sin AR) [sin (9-9. ) sin (pp, ) föß 
+ $ tan p, sin 4 Po ter] escl" 
where SB» = + sin" pa f 


The formula (4-21) is quite complicated, but can be 
simplified for systems not of world-wide extent. When the terms fof 
ana f58 are assumed to have a magnitude of 126,000, 000 or less, the 


terms containing the squares and products of (9-9. ) and GAM) may be 
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neglected in the factors by which the terms are multiplied. When 
these angles are such that the distance from the origin is 600 km., 
the squares and products amount to an order of magnitude of 
1:600,000,000 (33). 


The simplified formulas are:- 


(4-23) dE, {[sin (9-9. ) - 2 cos 9 sin p sin“ : (A-a )]8 B 


1 wc 
- 4 cos cos > (prp) sin > (9-9,)5 f 


3 ; e LE 
- (2 + i tan Po sin 49) sin (p 9,1182} escl 


(4-24) än. = [- cos 9 sin (X-A_)5 B 


l . "t 
+ y sin 9, sin n p, sin A) f5f] cscl 


4.3 The Use of Gravimetric Deflections of the Vertical. 


Probably the most accurate method of determining the 
equatorial radius from deflections of the vertical, involves a final 
correction by application of the gravimetrically determined 
deflections of the vertical ("absolute" deflections). 

The following analysis could have been given at the 
beginning of this chapter, but it was felt better to place it here, 
thus keeping the complete picture intact. 

Rice (25) p. 2, states that "a fundamental geodetic problem is 
to arrive at a reference ellipsoid of such dimensions and having such 


orientation that the deflections of the vertical are as small as 
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possible in the aggregate." 

An ideal reference ellipsoid would be centered at the center 
of mass of the earth and its polar axis would coincide with the mean 
axis of rotation of the earth. This is what Veis (32) p. 16, calls 
the Terrestrial Ellipsoid. <A further requirement which would be 
placed on an ideal reference ellipsoid would be that its surface best 
fit the surface of the geoid. 

The origin or initial point of a geodetic syster is 
generally assigned the astronomic latitude and longitude, p? and M^ 
observed at the place. The geoid distance N is generally considered 
zero and no errors in azimuth are assumed. The effect of making 
p= q^, and dj = A'at the initial point is a parallel displacement of 
the reference ellipsoid's system of coordinates with respect to those 


of the terrestrial ellipsoid and the geoid. At another point, P 


7? in 


the same geodetic system, which is displaced from the geoid, the 
normals do not coincide, the result not only of the displacement, but 
of the wrong computation surface, which is not the same size or shape 
as the true earth spheroid. The resulting meridian angle between q? 
and » is the astro-geodetic deflection (g’-p) in the direction of the 
meridian. (This discussion is generally in terms of this meridian 
component, for convenience, though similar remarks can be made 
concerning the prime vertical component, as well as the total 
deflection). 


The meridian component and the prime vertical component of the 


astro-geodetic deflection each consists of three angles in the same 
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plane, with respect to the normal to the mean or smoothed geoid. If 
we can correct this deflection in such a manner that both components 
are reduced to zero, what would remain would be the normal to the 
mean geoid. In doing so we have to change the size, shape, and 
position of the reference or computation ellipsoid so that our point 
Fi coincides (as nearly as possible) with the corresponding point on 


the mean geoid, and so the new ellipsoid is in the position of a 


terrestrial ellipsoid. The three angles are;- 


3 


(1) Correction for displacement, (a£, at.) 

(2) Correction for change of ellipsoid, (dE > dn.) 

(3) The absolute deflection of the vertical, or the angle 
between the true vertical and the normal to the 
mean geoid, or if preferred, the terrestrial 
ellipsoid of best fit. Since this is obtained 


gravimetrically, we call its two components z 


and io 
g 


See figures 3 and 4 for pictorial descriptions of this. 
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It should be noted that when the gravimetric deflections at 
the initial point are known, two of the unknowns of the previous 
problem are now known, namely ob s2 and ỒN + These are equivalent to 


the following: - 


(1520 ) A ~ geral, = rata 
Sans le na). - ea 
where (EL) x (o^ 3 Po) 


5. đ 
(Mad, = ($ - M) cos €, 


The gravimetrically computed quantities = and Mg? are in- 
dependent of the size of the reference ellipsoid (13) p. 239, even 
though they depend on gravity anomalies derived from the gravity 
formula employed. 

In the method (described in section 4.1) of determining the 
dimensions of the earth ellipsoid by means of correcting the astro- 


geodetic deflections with đến and đề 2 and subsequently reducing ES 


to a minimum, we have one effective method of solving the problem. 
However, we can go one step farther, as indicated in the preceding 
paragraphs, and make the final correction by subtracting the value Ea 
from the astro-geodetic deflections which have already been corrected 
by đến and dE. 

Thus our problem will be solved by a least squares solution in 


which 
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IS b minimum 


je 


Zune ng minimum 


In this way we are much closer to the mean geoid at the start. 
The result should be & terrestrial ellipsoid which best fits the 
surface of the geoid, providing, of course, that we have sufficient 
points from a large enough area for our computation. Theoretically 
the problem can be solved when astro-geodetic and gravimetric 
deflections are known at just two points, one at each end of a long 
arc, preferably as long as the earth’s radius. Naturally, the more 
points at which this data is known the better the solution. Heiskanen 
has suggested this method in a number of publications; see par- 


ticularly (12), and (13) p. 278-279. 
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Chapter 5. Flattening from Gravity Anomalies. 


A value for flattening can be gotten from the solution by the 
astro-gravimetric method described in section 4.3, but if this is done 
the gravity anomalies used in computation of the gravimetric deflect- 
ions are no longer correct since these anomalies are based on the 
adopted gravity formula, which in turn is based on a particular 
ellipsoid with its own flattening. This gives & clue as to how the 
flattening may be corrected from the anomalies themselves. see (13), 
p. 219. 

According to the theorem of Clairaut (13) p. 52, the 


following is true (extended form):- 


(5-1) 8= 272. 3n-2- iar 
where cM: are normal gravity at pole, equator, EP: 
m - (centrifugal force at equator)/? ,, = 7 
f = flattening Ñ 


The basic form of the gravity formula giving normal gravity 7 


for any latitude 9 on the ellipsoid is (neglecting the longitude term): 


(5-2) y * Yg (1 +8 since - € sin 2 0) 


where € = 2 mf - 5f 


e 
The problem is to correct the parameters Im B, and e, to the 
gravity formula where these parameters have been computed for a 


particular reference ellipsoid. Thus if B and c are corrected, a 
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correction to f is automatically found. The correction to e (based on 
currently used ellipsoids) is always so small, we may neglect its 


computation. 


Using, for example, the International Ellipsoid, the following 


values are used as & basis?- 


a - 6 318 388 m. f = 1/297 = 0.0033670 
y 978.049 gal m = 1/288.36 = 0.0034678 
B = 0.005288 € = 0.0000059 


* 


For obtaining a corrected value of 7, we have:- 


yayrthy 
where Ay x + y sin“ (p 
in which x' - unknown correction to Vp in mgals. 
y’ = unknown correction to B, in mgals. 


Before correction, the gravity anomalies, based on the 


ellipsoid used and isostatically reduced, are:- 


PB - 8S - y 
After correction they are:- 
= = a= = 
VIENE y c MM eo) 
6-7) =- A7 R A 


Then, we can write error equations for each point where 


gravity observations have been made:- 
‘ pee 
x” + ysin @ = 4=-v 


If we have enough such stations, spread Out over a large 
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enough area, we can solve the problem effectively by the method of 
least squares in such a manner that:- 
2 seme 
25 Vv = minimum. 


An example of this type of solution made by Heiskanen, in- 
volved using anomalies from over 1500 squares of 19 by 1? àimensions 


(13) p. 77. The resulting parameters were:- 


In = 978.0451 
B = 0.0053026 
€ = 0.0000059 


The flattening corresponding to the corrected B was:- 


f = 1/298.2. 





Chapter 6. Gravimetric Deflections of the Vertical 


6.1 General. 

The gravimetric method, of which the subject of Cnapter 5is a 
part, also produces what are sometimes referred to as "absolute" 
deflections of the vertical, obtained from gravity anomalies. The 
absolute or gravimetric deflection of the vertical at a point is 
defined as the angle between the normal to the geoid and the normal to 
the earth spheroid at the point. -- Veis (32) uses the term "mean 
geoid" to describe the surface of the earth spheroid in the area of a 
point, meaning a smooth surface striking a "mean" between undulations 
of the geoid. -- The discussion in Chapter 4 (section 4.3), describes 
the use of these deflections in determining the size of the earth 
ellipsoid. This chapter will present a brief outline of the method of 
determining the deflections from gravity anomalies. 

To begin with, the gravity anomalies Ag, are caused by dis- 
turbing mass anomalies, and with this in mind and with the aid of 


fig. 5, we can easily see how local absolute deflections occur. 


Geoid 





< 
NI Sphero¡d 


Normal # ; 
Du S Vertical 
Io) 
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The mass surplus tends to "pull" the plumb line toward it, and 
conversely the mass deficiency tends to repel it, thus causing the 
deflection £ = P' - Y. The sign of this deflection is apparent from 
this expression which shows that £ is always referred to the 
astronomic latitude, q', of the point in question. The sign conven- 
tion for gravimetric deflections is the same as for astro-geodetic 
deflections. Referring to figure 5, we wee that when 9 > g', £ = 
(p' - @) < 0, and when 9 < ọ', £ = (m' - p) > 0. In U.S. practice A 
is considered positive westward, and in European practice positive 
eastward. In either case, when X>A', N= (A' = X) cos M < O, anà 
when X <A", P = (A! - A) cos p > 0. 

From a theorem developed by Stokes, Vening Meinesz formulated 
expressions for absolute deflections of the vertical in 1928.  Stokes' 


Formula expressed the geoid distances N, as follows:- 





eT TT 
(6-1) v - E fa a | sw sin y Æ d y 
O 


O 


in which y = angular distance of a circular ring of width dy. See 


Pee 6, 
Ag = average gravity anomaly in ring, having the limits y and 
y + dy. 
R = Mean radius of the earth. 


G - Mean gravity. 


f(y) = 5 S(y). 


(6-2) S(y) = ese > Y + 1-6 sin > Y - 5 cos y - 
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- 3 cos y 1n (sin 5 Y vaain" L y) 


which is the original Stokes Function. 


sin y dy QA = surface area element for unit sphere. 


POLE 


P= Computa tion 
Point 


C = Center of 
Sphere 





AS = Surface Area 
Element 


Fig. 6 


Vening Meinesz got his formulas by differentiating the Stokes 


formula, as indicated in fig. 7, in the form:- 


is 
(6-3) N= 57 OR Jr) ages 


ORT 
Ba 77, ^N Af 


£ 
dx 
N 
T—— d^ 
Fig. 7 | 


Upon differentiating, the following are obtained: - 


SPHEROID — 
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"t 
(6-4) dp" E Jem cos A^g d S 
2TGR 
ft 
(6-5) ann = E (y) sỉin A Ag đ 6S 
euch 
where dS = R^sin y d y dA, and 


1 
cos Ir 
(6-6)f'(y) = zl = =- ——— + 8 sin y - Ó cos 2 
2 sin 2Ý 


o 


- 3(1 - sin 2 csc y + 3 sin y 1n [sin 2 (1+ sin sv)] 


Application of these formulas to practical use requires 
further manipulation. They are used in different form, whether one is 
going to compute the effect of areas beyond y x 3°, or within the area 
circumscribed by a circle of radius y = 30 (39 is an arbitrary dividing 
line, which, though close to the limit, has been chosen for convenience 
(31) p. 96). These different forms are described as the "Circle -Ring" 
Method, for use up to y ~ 3°, and the "Squares" Method, for use 
beyond. Still another form is used for the small inner circle close to 
the computation point, having a radius ro generally, of 5 km. or less. 
Figure 8 shows this division of methods pictorially. For further 


references, see (13) p. 63-70, 257-264, (15). 


6.2 The Squares Method. 

This method, used for the earth represented as a spherical 
surface, and for the area beyond a radius of about 39 from the comput- 
ation point, makes use of mean gravity anomaly values for each P by 
T square (spherical trapezoid) between y = 3° and y = 20°, and for 


each 5° by 5° square beyond 20° and around to the antipode (31). 


h2 








5x5 SQUARES 






= /°x/° SQUARES 


ZZ 






I = INNER CIRCLE 


Fig. 8 


These should be either free-air or isostatic anomalies, the latter 
being best. If Bouguer anomalies are used they must be corrected to 
obtain free-air anomalies. To adapt formula (6-4) to the squares 
method, we first replace dS by R“sin vavyd A, obtaining:- 


esc 1" 


(27) de" = E (y) cos A^g sin d vd A 


Then using the square element dq on the unit sphere, R = 1, 


(14), we obtain:- 





(6-8) dae" = DL few cos AAg dg 
q 


The expression for dn" being the same with sin A replacing 
COS A. 
When Ag is in milligals, and when G is taken as 979.8 cm. /seez, 


the term before the integral sign will be: ‹ 
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- 070335 
One procedure is to perform numerical integration summing up 


the effect of each square over the earth as indicated in the following 


formulas - 
m "t E tt 8 
(6-9) dE" WT CR cos A, 28, 4 


where q is the area of a da by 1? or 59 by 5° square and depends upon 
its mean latitude. The subscript q, where occurring, signifies mean 
values of the quantity for the square in question. 

The squares method adapts readily to solution by high speed 
electronic computing machine, but manual graphical solution, though 
more time consuming, is made easier by use of maps on polar stereo- 
graphic projection. With this projection templates can be easily 
prepared for use on the large area, small scale maps required, since a 
feature of the projection is that all circles on the sphere project as 
circles on the map, so azimuth lines (which are great circles) and the 
circular rings (small circles) can be easily drawn. Necessary 
formulas for preparing such templates may be found in (13) p. 259-260, 
and in (31) p. 98-100. 

The functions f'(y), f'(y) sin v, `... sin y d y, have 
been computed by Sollins at close intervals for different angular 
distances y, from the computation point. These tables have been 
published by Sollins (29) and are very useful in the solution of the 


problem. 





Tn 
6.3 The Circle-Ring Method. 


In the relatively close &re& around the computation point 
(within y e 3?), where a detailed gravity anomaly map is available, 
formula (6-7) is used in the following form:- 
y+dy 
„n ese 1" , 
dt er cos d af £* (y) sin y ^8, d y 


Ay y 


(6-10) 


y+ay 
esc 1" 


= 2TG (sin A, - Sin e P (y) sin y ^8 à y 


A 


For d n", substitute (cos A, - cos Ay) for (sin A, - sin A). 


Formula (6-10) gives the effect of a circle-ring compartment bounded 


by y and y + dy and by A and Aus and having a mean anomaly Ag, , on 


Y 
the total component e 


T 
Values or / f'(y) sin y d y are tabulated in Sollins Tables. 


4=10m y+dy 
Thus for the value of a particular / £*(y) sin y à y, the tables 
y 
would be used as indicated in the following expression:- 
y+dy 
| £* (y) sin y d y - few sin y d y. 
y=10m y= 10m 


Figure 9 gives pictorial explanation of the components of the 
method and corresponding templates. 

Any system of azimuth sector and ring zone intervals can be 
used in this method, but a system which yields equal effect for each 
compartment would obviate much unnecessary work. Kasansky developed 


a system in which the sine differences are 1/8 (32 compartments), and 
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COMPARTMENT 





P = computation point 


y = inner radius of compartment 
y + dy = outer radius of compartment 
A, > A, = bounding azimuth lines 
A Ey = mean gravity anomaly in mgals, read from 
iso-anomaly lines on map 
Figo 9. Components of Template for Circle- 


Ring Method 





the ring radii, beginning with Da 5.0 km., follow the geometrical 
progression = 1.270 r_. Farther from the station, the ratio is 
slightly smaller than 1.270 because of the earth?s curvature. Tables 
of these Kasansky Circle-Rings may be found in (13) p. 267, and in 
(16) p. 17. The effect of each compartment is 0.001 radial deflection 
effect per 1 mgal Ag. 

D. A. Rice, of U.S.C.& G.9., published a system of Circle- 
Rings (25), in which each ring was divided into 36 equal azimuth 
sectors of 10^ aperture. Heré also, the radial effect for each com- 
partment for 1 mgal Ag is 0°001. Rice’s method for computation of the 
radii involved using the formula for the deflection effect along any 


radial lines- 


ese 1" 


15-11) e" = Ds Feo sinvdvZgdA 





so de” = wo" cos A 
an" = d 6" sin A, and 
ese 1" Hạ 
(6-12) de" = i E ty AA 
Yn — 
where ly is obtained from Sollins? Tables (29). 


Then for AA = 10° = 0.17452 radians, and for G = 981,000 mgals, 


formula (6-11) becomes: - 


y 
ae" = 0.0058405484 tì Age 


O 
For d0" = 0,001 and Ag = I mgal, 


AT. = I = 0.17121680 


a “ a ¿01 i “eet 
. re ee 
Fe ee Ì 


— - a $ 





giving Al, the circle-ring table was prepared. 


Vo 1 
I -I. = 0.17121680, etc. 
ie Yo 


Then selecting interpolated values of y from Sollins' table 


Table I is Rice's 


table reproduced. 


Table I 


Rice Template Radii for Gravimetric Deflection of the Vertical 


E 


bere 
& I D o0 0-1 OW £u Io P. 


For radial deflection effect of O7001.and mean com- 
partment anomaly of l mgal; angular aperture 10. 
(n » zone number; r = inner radius in km.) 














r um r A. r = r 
0.119 14 1.099 27 10:15 39 T-T 
0.141 T5 1. 30h 28 12.05 LO 92.22 
0.167 16 1.547 29 1^.29 4] 109.0 
0.198 17 1.836 30 16.94 42 128. 
O.235 18 2.179 31 20.09 43 151.9 
0.279 19 2.586 32 23.83 ul, 179.1 
0.331 20 3.068 33 28.25 45 210.9 
0.393 21 3.641 34 33.48 46 248.0 
0.467 22 4.320 39 39.67 47 291.2 
0.554 23 5.125 36 47.00 48 341.2 
0.657 24 6.081 3T 55.66 ho 399.0 
0.780 25 7.216 38 65.90 50 465.5 
0.926 26 8.560 


With the use of circle rings such as Rice's, the problem of 
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computing the gravimetric deflections becomes simple, though laborious. 


All that is necessary is to estimate the mean anomaly from the contours 


in each compartment, make summation of Ag for each azimuth sector, mul- 


tiply by 0.001 cos An or 0.001 sin A? and the effect for that sector is 


solved. 


effect. This may be further simplified when one takes 


note of 


Then the effects of all 36 sectors are added together for total 
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the fact that mean azimuths in 180° opposite sectors will have the 
same cosine or sine, but with opposite sign. One must be careful with 
signs, since both the azimuth of the sector and the anomaly gradient 
contribute to the sign. Additional information on this aspect of the 
problem will be given in the explanation accompanying the computations 
made in Chapter 7. 

In the event gravity material is poor, which is too often the 
case (and one has to be a good judge concerning when poor material is 
to be considered good), Heiskanen (13) p. 268, recommends combining 
two zones and two sectors into one four-compartment compartment, the 


effect of which then would be 0004 per mgal. 
6.4 The Inner Circle. 


The circle-ring method, although it can be applied beyond 
y = 3 should not be used very close to the computation point since 
the term csc 5 Y in the value of f’(Y) sin y, in formula (6-10) 
approaches infinity. Rice (25), suggests a lower limit of 100 meters. 
However this is usually too small because of the sparsity of gravity 
anomaly values, except in isolated instances, within any given 100 
meter circle. 

The inner circle, being small, is considered to have a 
uniform anomaly gradient, so the anomalies read from the cardinal 
points are used in tbe following expressions, the derivation of which 


may be found in (31) p. 93, 100-101:- 


-Ầ.............. a 6 (eg) x 
(6-13) de" = ET pe zur, = 0.105 ro 5 x 
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an" = SCH xà Y y = 0/105 r_ 2 (48d 
o 6y 


where E radius of inner circle in km. 
5 (48), = 28, - 48, in mgals (N. hemis. ) 
5 (48), = 02, - 4g, in mgals (European practice) 
8X ST OV 222r 
O O 


Substituting these where appropriate, formulas (6-13) become: - 


(6-14) ag" = 0/0525 (48, ”- ^g) 
an" = 020525 (48, - 48,,) 


See figure 10 for further clarification of this "one- 


gradient" method of computing the effect of the inner circle. 


/ 


Ag. 





Fig. lO. Inner Circle - One-Gradient Method 


Rice used & three-gradient method rather than one-gradient as 
above, in computations for 16 points in the U.S. (26). Taking 


anomalies at NE, SE, SW, and NW points, as well as at the cardinal 
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points of the circle, three N-S gradients and three E-W gradients were 
used, with weights of one-half given to the outer gradients and a 
weight of one to the center gradient. Figure 1l illustrates this. 
This method should be used only when detailed gravity information is 


available inside the inner circle. 


Ôx (S-N) 
dy (W-E) 


Ox (SW-Nw) 
(SE-NE) 

by (Sw-5£) 
(NW-NE) 


/. 4144. 





Fig. ll. Inner Circle - Three-Gradient Method 


" 8, - “En 
de"(S-N) = 0°105 y I. Weight = 1 
Man 7 48 
ag" (SE-NE) = 01105 5 r Weight = 5 
Bea, “28 
ae" (Sum) = ov105 —" r , Weight = 5 
X Weights = 2 


Combining these, using weights shown and values of ồx given 


in fig. 11, we obtain the following expression, given in (31) p. 101:- 


(6-15) ag" = 0.0262 (28.- 287) + 0.0186 (08 .-08,.+08..,-28,5,) 
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This is for the northern hemisphere. For the southern 
hemisphere, just reverse S and N where occurring. The other component 


of the deflection is similar:- 
(6236) a"! = 0.0262 (02, ,-087) + 0.0186 (^g... "€. s PE ^w) 


This is for European practice, in which X and n are considered 
positive eastward. Where the practice is opposite, W and E should be 
reversed where occurring. 

As may be deduced from the opening paragraph, the smaller the 
inner circle, down to certain limits, the more concentrated the gravity 
information must be. Vening Meinesz (13) p. 70, states that the error 
incurred using a circle with r, - 30 km. is less than 1%. A good 
round number seems to be 5 km., and for good results there should be 


from four to eight points of known Ag within the circle (13) p. 267. 





Chapter /. Computations 


7.1 Dimensions from the 98th Meridian in the United States. 


(From Astro-Geodetic Deflections ). 


Using data from U.S.C. & G.S. Special Publication No. 229, (TY 
consisting of astronomic and geodetic latitudes and meridian com- 
ponents of the astro-geodetic deflections of the vertical for 30 
points between 97° W and 99° W longitude, extending from North Dakota 
to Texas, values for a and f will be determined. The astro-geodetic 
deflections as published in the publication named have been reduced 
according to the Hayford isostatic assumption for a depth of com- 
pensation of 113.7 km. 

The problem to be solved is E = = minimum, according to the 
theory of least squares. The Vening Meinesz displacement and change of 
ellipsoid formulas will be used. To adapt them to practical use, terms 
containing the same unknowns must be collected and combined. 

Prior to solution, a test problem was solved to determine what 
terms may be disregarded and whether or not the long or short formulas 
should be used. The maximum difference in longitude is about 1°, ana 
terms containing longitude were determined to make negligible con- 
tribution to the residual deflection. Also, for the arcs of limited 
length (AQ = 1373, in this case) encountered in this problem, the 
difference between residuals computed by long and short formulas, 


turns out to be negligible. Terms containing 5N will also be 
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neglected to simplify the problem. 
The Vening Meinesz formulas, given as (4-19) ana (4-23) in 


section 4.24, are here repeated with terms omitted as noted above. 


Displacement: - dé, = cos (o, > Po) o6. 


Change of Ellipsoid:- dé p,58 - p.df - pif df 
c 1 2 5 


where 8B = 5a sin“ 9. ôf 


Substituting the foregoing expression for 58 in d£, and collecting 


a 


terms, we obtain: - 


da e 
dE. = Py et p, sin p of - pof = pyiof 
= A + [ i E £j Of 
a Pay. ea E 
in which p, = sin (o, - 9.) 
p, = 4 cos M cos (o +p ) sin Io - Q) 
2 n UNTEN ern O 
2 3 
Pz = (2 + y, tanp_sin ho) sin (o, - 9,) 


The observation equations will take the form given below, with sub- 
script O signifying the initial point, Meades Ranch, Kansas, and sub- 


script n signifying any other point in the arc. 


EL cea E (AO 


For convenience we transform the unknowns and the coefficients 


as indicated here: - 


ure 
1 


5 a 
(h) dg. + (3) 1000 + (k) 10,000 8f + (1) 


hx + jy + kz +1 





where 


and 


1 
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OE 


O 
HN unknowns 
1000 y 
10,000 Of 
cos (9, - 9.) 

à 1000 esce 1" 
sin (9, - 09) 85 — 


| sin (ọ -0 )[ sản” @ - (2+È tan 9. sin k 9 )f] 


1 n. cese 1” 
- ^ cos.Q, cos 5 (p +P) pud (9-9) ] 10,000 





Table II - Basic Data - 98th Meridian 


Initial Point:- 


Meades Ranch 


— 
— 


1/294.98 
0. 00339006 


v. - 39713'067686 N. 


Point Latitude Astro-Geodetic 
n Location Ph Deflection 
(p” - q) 
406 N. Dak. 46%02'17"N. Tag 
571 S. Dak. 45 28 16 1.06 
566 t 44 53 53 0.78 
569 i hh 00 52 -2.93 
570 ý 43 43 24 -1.43 
567 a 43 18 48 0.96 
412 Nebr. h2 25 29 3.88 
407 " 41 35 44 slk 29 
413 : ho 46 35 -2.93 
235 Kans. 39 h6 35 -2.03 
Meades ü 39 13 2T 0.00 
238 i 38 43 47 -1.29 
245 i 37 51 38 -2.51 
240 P 37 14 28 20910 
522 Okla. 36 ko 08 0.20 
524 i 35 56 49 1.88 
516 u 3586 27 ge T). 
521 : 34 38 53 1.98 
602 Texas 3937 21 -1.28 
620 i 32 48 10 -0.91 
639 p 3215 45 6.78 
605 i 31 39 26 2.20 
594 i 30 59 h6 (665 
592 , 30 16 21 =), 
635 y 29 42 49 -4,.05 
652 í 29 22 16 2.14 
627 $ 28 52 38 -0.59 
634 \ 28 18 06 -0.68 
288 : 27 44 35 1.44 
630 : 27 ho 16 4.91 
Clarke’s 1866 Ellipsoid:- a = 6,378,206 m. 





Table III - Observation Equations - 98th Meridian 


Point h j k 
n (x) (y) (z) - (9. -9,.) 
406 0.993 3.837 ale Se a 
571 . 994 3.519 232850 3806 
566 . 995 3.19T -1. 348 0. dfe. 
569 997 ROL -1.17 -2.93 
570 4997 2,995 «128193 al. 43 
567 - 997 2.306 - 14088 0.96 
412 .998 » 1.806 -0.835 3.88 
407 . 999 1.338 = 3687 -1.29 
413 1.000 0.876 - .429 -2.93 
235 1.000 .312 e WAS -2.03 
Meades Ranch 1.000 -- -- 0.00 
238 1.000 - .279 .146 -1.29 
245 1.000 - „770 .414 -2.51 
240 . 999 -1.119 .613 -0.40 
522 . 999 -1.442 . 903 0.20 
524 . 998 -1.849 1.058 1.88 
516 . 998 -2.228 1.290 dir 
581 . 997 20530 14520 1.98 
602 . Q95 -3.157 1.911 21239 
620 . 994 -3.617 2.367 -0. 91 
639 .93 -3.920 2.457 6.78 
605 . 991 -4,258 2.710 2.20 
594 . 990 -4.628 2.993 7.65 
592 . 988 -5.032 3.310 AL 
635 . 986 28898 30567 -4.05 
652 . 985 5. Qe 3.716 adil 
627 . 984 -5.808 3.944 -0.59 
634 . 982 -0 129 4.214 -0.68 
588 . 980 -6.437 4.481 1.44 


630 .980 -6. 4TT 4.516 4.91 





Table IV - 


Forward Solution:- 


x y Z 
x 429.62 -17.43 +35. 76 
-1 +1.601263 -1.207292 
y +398.56 -250.22 
- 75.95 + 57.26 
(1) +322.61 -192.9% 
-l +0. 598122 
Z +160. 48 
- 43. 17 
-115.41 
"m + 1.90 
- 1 
1 
[££] 
Back Solution: - 
Z Check y 
il -1.131579 -2.131579 +0. 199622 
Z -0.676822 
y 
> -1.131579 -2.131579 -0. 477200 
E(check-U) -1.000000 


Standard Error of single observation?- 


_ M LEE], rá 200.02 = + 


least Squares Solution - 


1 > 
+15 08 -33.02 
-0.509115 -1.115125 
88.55 412.36 
425.15 452.89 
-64.40 465.25 
+0.199622 -0.202257 
+58, 88 + 4.90 
«18.21 -39.88 
-38.52 +39.03 
+ 2.15 + 4.05 
-1.131579 -2.131579 
4223.62  +209.03 
- 760 <= 107592 
- 12.86 + 13.03 
= 2 ...... 
+200.65 4200.66 
Check X 
-0.202257 -0.509115 
-1.274944 +1. 366146 
-0. 764132 
-1.47720) +0.092899 
-1.000001 
7.43148 =+2!'726 


S 


98th Meridian 


Check 


-1.115125 
+2.573438 
-2.365417 


-0. 907104 
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Table V - Computation of Correlate Numbers - 98th Meridian 
1(x) l(y) E) 
(x) + 1 
(y) O 
+1.6013 
+1.6013 TT 
(z) O O 
-1.2073 +0 
+0. 9578 +0.5981 
-0.2495 +0.5981 + 1 
1.90 O = HỆ 1.90 Qe = 0.5981 1.90 e = -0.2495 
Q,, = 0.5263 qj; 7 0.3189 Q, 7 -0.1313 


XZ 


322.61 Bes" 192.96 2 = 1 
322.61 Ay = 1 + 60.7438 = 61.7438 


Q y = 0.1914 

322.61 A - 192.96 x” 1.6013 

322.61 x = 1.6013 - 25.3356 = - 23.7343 
Q x - 0.07357 

29.62 Qe - 47.43 os + 35.76 om 

29.62 oms l - 3.4894 + 4.6953 = 2.2059 
E = 0.07447 


Standard Errors of Unknowns after Adjustment: - 

2.726 Y 0.Ó7447 = 2.726 (+ 0.2729) = + 0.7439 
2.726 J 0.1914 = 2.726 (+ 0.4375) = + 1.1926 
2.726 J 0.5263 


il 


3 
u 
T 


AE 


2.726 (+ 0.2294) = + 0.6253 
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Results of least squares solution:- 


x = + 0.0989 = ot 
n _ ba 
y = - O.4}TT200 = 1006 
Z = - 1.131579 = 10,000 8f 
EE - NOSOP OA 
Sa = 1000 (- 0.4772) = - 477.2 m. 
BET 
Sf = —316,000 4 ^ 0. 000113158 
a, = 6,378,206 m. 
da = -477 
New a = 6,377,729 m. + 1190 m. 
f,® 0.003390060 
ôf = -0.000113158 
New f = 0.003276902 + 0.00006253 
1/f = 304.9 + 5.4 


The solution, as can be readily seen, does not conform to what 
is known to be the truth, and the standard errors resulting show that 
the precision of the solution is "somewhat" rough. This of course can 
be explained by noting that the length of the arc is quite short, and 
that the number of points used is small. Better results certainly 
could have been obtained if many points having greater longitude 
differences had been used, thus affording the benefit of the prime 
vertical components of the deflections. However, the method used is 


what is intended to be demonstrated here. 
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It could have been expected that the results obtained from the 
data used would be as they are. Differences in latitude between the 
initial point and other points in the arc apparently differ systemat- 
ically. It is seen as a rough check of the standard error of the 
value of a obtained in the least squares solution, that if we take the 
average value of (9’-~) for the five northernmost stations and of the 
five southernmost stations, the resulting difference between the two 
is - 8.57. For the arc as short as the one used this is high. 

Taking the mean latitudes of each of the two groups of five 
points, we get Ap = 16°26’. Taking the length of A (q’-) at 31 
meters per second of arc, we have - 8:57 @ 3l m. = -266 meters. Using 


one radian as 9985. we can determine the proportional difference for 


an arc of the length of the earth’s radius:- 


uu ( - 266) = - 930 meters. 


o 


This kind of rough calculation could tell us in advance 


approximately what results can be expected. 
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7.2 Dimensions from Sixteen Stations in South Central United States 


by the Astro-Gravimetric Method. 


For an illustration of the astro-gravimetric method as 
described in section 4.3, the published data on 16 points in South 
Central U.S. (including Kansas, Missouri, Oklahoma, Arkansas, Texas, 
and Louisiana) as computed by Rice (26), was used. 

In the referenced publication, all geodetic data were con- 
verted from Clar’e’: 1866 Ellipsoid to the International Ellipsoid. In 
addition, as a result of the investigations, Rice determined that 
several of the astronomic observations which had been published in (7) 
were wrong, and were consequently done over. Moreover, all latitude 
and longitude determinations were reduced to the ss system of 
cataloging mean places of stars. The gravity anomalies used were free- 
air. 

All the data needed for solution of the present problem was 
given in Tables VI and VII of ref. (26), namely the latitudes and 
longitudes of each station, as well as the astro-geodetic and the 
gravimetric deflections for each. Here again, as in the problem of 
section 7.1, the extent of the area covered, and the number of stations 
involved are inadequate for a meaningful solution to the problem, but 
for illustrative purposes, the results are gratifying. 

Meade’s Ranch, Kansas, was again chosen as the initial point. 
The Vening Meinesz abbreviated formulas were determined to be adequate 
and were used for both £ and "n. Since Meade's Ranch was one of the 16 


points given, the values of OEA and on are known, thus reducing the 
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number of unknowns to two (EN, again being neglected). 
In adapting the formulas (4-19), (4-23), (4-20), and (4-24) to 
this problem, the terms were grouped in the logical manner for ease in 


computation of the coefficients as follow:- 


Ca M 


E = 8/5, +h,09, +3, ux +k, (10,000) 8f * (9 -9) - En 
g} = + cos (P - 9) : 
h, = + sin q sin GÀ) 
100 
je [sin (p - 9. ) - ? cos 9 sin 9 sin’ 5 (A EC )] ne 1" 
j ] 2al 2 
k =+ | (sin (o - 9. ) - 2 cos 9, sin q sin 5 2] sin Q. 


1 
(+ p sin 5 (9 - 9,) 


Mle 


- 4 cos @ cos 


3 esc 1" 
- (2 + T tan Po Sin + Po) sin (p - a ioo 


n = an, + an, + (A' = dX) cos Q -N 2 


TS 8,8 Ra ve h, 8n. ds = + ko (10,000) 5£ + (A' - A) cos Q - ng 


g = - sing, sin (A - ho) 
A, = 4] 

100 " 
jo 7 - cos 9, sin (A- A o) me SEM 


k, = [ - cos Po sin (A - No) sin 95 


csc 1" 


1 : : 
4 n sin ps sin 4 PE sin (A- No? 2 10,000 





dị T65 


n = 
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Then the observation equations become: - 


tk (10,000) 8f + 1 


1 
where 1, = 8, 3, + hd Dept en = - 
E79 -0 
S to "uelis - See `”. 
6n cNEM - Bo = a - 
D = + k, (10,000) öf + 1, 
where 1, = 8,85, + hồ "+ Ma = n.) 
n = (A' A) cos 9 


a 
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Table VI - Basic Data - Rice's 16 Points 
Point dime p X 
n 
i Meades R. 39 13" 27” y. ‘G8 32" a1" xw 
2 Arcadia 37 38 08 gie 35 X 
3 Dirks 37 25 HỘ 99 16 21 
4 Twin 36 06 57 96 46 ho 
5 Little Rock 34 4h 51 oo 16 25 
6 Bogue 3 3 05 90 23 1H 
7 Ecore 33 36 20 92 54 25 
8 Polk >> 23 O5 99 05 16 
9 Bartley " EB 12 0 101 ho 44 
10 Roby 32 4k 36 100 22 46 
11 Lacasa 32 39 05 98 u 38 
12 Burns 32 35 28 93 14 20 
13 Sears E? 33 32 100 02 16 
14 Brooks 32 25 09 Je OT “Ol 
15 Bynum 32 19 16 lol 00 00 
16 Legion 32 16 19 90 06 07 


International Ellipsoid:- a = 6,378,388 m. 
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Table VII - Observation Equations - Rice's 16 Points 


Point j k 1 

t 1 0 O O 
2 - 0.093 + 0.672 + .3 
3 = Ace + @773 = SD 
4 "M + 1.366 4» M6 
5 - 261 + 1.993 + 1.54 
6 < 5223 + 2.484 + 1.44 
7 - 48293 + 2.554 + 95 
8 - 329 + 2.681 T 5 
9 = "¿UL + 2.767 + 02 
10 - $360 * 3.001 mm 55 
l1 -. OO + 3.050 7 .o4 
12 =. a379 + 3.065 + 1.44 
13 = ame + 3.0% + 39 
14 - .392 + 3.146 + .8 
15 - .390 + 3.215 + .45 
16 - 406 + 3.205 + 1.31 

N 1 O O O 
2 - 0.172 + 0.440 - 0.90 
3 + .032 eu S 20205 
4 =- .O77 + 197 - .29 
5 - 275 + .702 + .36 
6 - 334 + .350 + .33 
" - ab 1 627 + 1.12 
8 + 024 - ,@62 + yer 
9 + «137 - .3h9 + 2,53 
10 + .080 = 2205 + 1.02 
11 + 007 = 3017 + 1.15 
12 - 288 + .590 + .60 
13 4+ 005 - 167 + 1.43 
14 = 28C + .714 + .83 
15 + SO aa + 2,16 
16 335 + .937 + .65 
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Table VIII Least Saueres Solution ~-  nice's 16 Points 

BR 10,000 8f 1l 2 1 (a) 1(f) 
100 : 

òa_ 

100 + 2.17 - OR - EOS =u15266 l 

- 1 +6.54839 | «1.66820 +7.21659 
10,000 8f + 106.07 + 927.95 + 119.81 0 
- 93.09 - 2207) 4302959 09549 


13.02 + 4.2: + 17.26 +6.548 1 


+ 


- 1 -0. 3256510 1.32565 

T > + 30084 + 55.17 

- 6804 => "BEN? 

1026 A FRO? 

[vv] + 03,42 7+ $3.53 
da 


10,000 Of Check Check 


100 


l  -0.32565 -1.32565 +1.66820 +7.21639 
10,000 èf -2.13248 -8.68087 


ES 7?20:39565 -1.32505  -0.4042O9 "1, HOURS 


Standard Error of Single Observation:- 


_ AI Peete WW. 
EN Me Xu = (0889 
"...”. 
mes ^ 099 
(+ 2.17 Oia PERS TE 
O a 
Maa ” o c cm 
4 ES 
Br uo ^ 999 
Standard Errors of Unknowns;- m= pN Q 


m 
= 0.8836 Y 3.755 = 1.712, 10,000 m,=0.8836V0. 07 =0.24520 


m = 171.2 m Mp = 0.000024520 
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Results of Least Squares Solution:- 


10,000 8 f - - 0.32565 + 0.24520 
ö f = - 0.000032565 + 0.00002520 
£, = 0.003367003 
f = f +8 f = 0.003334438 + 0.000024520 
== 299.9 + 2.2 
5 a 
100 = - 0.46428 + 1.71 meters 
Sa = = hóa a 
a, = 6,378,388 
aca tba = 6,378,342 + 171 meters 


Here, as in the problem of section 7.1, the standard errors 
of the corrections are almost as large as or larger than the 
corrections themselves. But the fact that they, as well as the 
corrections to f and a themselves, are relatively small, testifies as 
to the value of the method. However, it must be realized that the 
standard errors reflect only the precision of the least squares 
adjustment, assuming the given data has no error, which of course is 


not true. 
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7.3 Computation of the Flattening from Gravity Anomalies. 


Given a set of isostatic gravity anomalies (Hayford, 
D = 113.7 km.) at 67 points at various random positions scattered 
over the earth, and computed on the basis of Helmert's gravity 
formula of 1901, we can demonstrate the solution of a problem of 
correcting the gravity formula and of obtaining a value for the 
flattening of the ellipsoid (13), p. 76-79, (14). The method is that 
of Chapter 5, but extended to.include the longitude term of the 


formula. 


Helmert's Formula; - 


2 


Y = 978.030 (1 + 0.005302 sin q + 0.000007 sin“2o) 


for which the B term yields the value of f - 1/298.2. See 
(13), D. 72. 
The general form of the gravity formula, including the 


longitude term is:- 


y Y Ll +B sino +E sin 2ọ + r o cos 2(X - A )] 


Yall + B sin presin 2p+r cos^o (cos2Acos2X + sin 2X sin 2x. )] 


Making corrections, we have Y! = 7 +A) , 


where Ay = x' + y'sin^o + z!cos^ cos 24 + u'cos^g sin 24 


x' = unknown corr'n to Ye in mgals 
E N i Ñ B in mgals (+ 1.022 x 1073) 
NE y N " “Fees 24, in mgals (+ 1.022 x 10°) 


Du á 2 +... 24, in mgals ($ 1.022 x 1073) 





(0 


Using /Y - /8 = - Vy we have 
: - : 2 ; 2 : 
Xx" + y'sin q + 2'c06 0 cos 2 A +m'cos 0 sim 2 A ¬-‹@@ = - V 
The observation equations will then be:z- 


ax'+by'+cz'!+dAu!-1=-v 


where a = 1 
b = sin p 
2 
c = cos "@ cos 2 X 
d = cos“ @ sin 2 À 
1 = & 


The solution, according to the method of least squares will 


be E = minimum. 





Table IX -  Isostatic Gravity Anomalies (Hayford, D - 113.7 km.) 


Point p A ^g Point p A Ag 
1 1995 1195 +16 34 2295 3695 4 o6 
2 79.5 17.5 + 8 35 JS 38.5 + 48 
3 T9 15.5 =- 26 36 14.5 - 75 +44 
4 71.5 DODENUS 37 13.5 42.5 + h3 
5 65.5 24. 5 38 - 1585 38.5 + 27 
6 64.5 liege 28 39 - 2685 15.5 + 34 
fi 62.5 17.5 -20 Lo po) 1260.6 >- 1 
8 60.5 22.5 + 2h 41 49.5 94.5 + 9 
9 58.5 38.5 +16 42 48,5 116.5 +10 

10 67.5 40.5 + 28 43 47.45 06.5 +19 
11 55.5 13.5 + 10 uh 47.5 795 + 6 
12 4.5 39.5 E 11 45 46.5 123.6 - 5 
13 54.5 55.0 ıı 46 46.5 92.5 +58 
14 53.5 290.5 +15 47 46.5 6645 >- 1 
15 93.5 50.5 + 23 48 45.5 102.5 + 23 
16 5275 25.5“ 15 49 44.5 110.5 + 14 
17 5 - 0.59 9 50 44.5 714.5 + 26 
18 9.5 11.5 +36 51 43.5 ob. p + I4 
19 Dor 10.5 32 52 43.5 76 - 16 
20 48.5 16.5 + 31 53 4275 71.5 + 1 
21 7.5 39.5 +13 54 39.5 1205  - 20 
22 16.5 TIE DO 55 39.5 ¿LN - 1Ì 
23 16.5 34.5 +6 56 38.5 7949 - 5 
oh 45.5 10.5 + 30 5m 3695 87.5 + 1h 
25 HiS- = 1.55 - 58 355 092 +2 
26 44,5 39.5 +31 59 34.5 79.5 +10 
27 43.5 12,5 +13 60 315 94.5 =- 4 
28 40.5 14,5 + Th 61 56.5 59.5 +55 
29 39.5 53.5 +15 62 5145 104% - 2 
30 36.5 + 6.5 +14 63 46.5 61.5 + 1h 
pal 3545 14.5 + 70 64 38.5 76-5 - 66 
32 28.5 33.5 + 17 65 23.5 79.5 + 18 
33 26.5 33.9 + 31 66 22.5 mo + 

67 line 78.5 - © 





Table X 


Point 


W 0-3 ODM FU NN P 
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- Error Equations - Flattening and Correction 


PPPRPPPPPPPPPPPPPPPPPPPPPPPPPPPPrPrPp 


b C 
y) (z' ) 
. 9669 . 0306 
. 9669 . 272 
. 9456 .0467 
.8993 Si 
.8281 > .1128 
. UA . 1706 
. 7868 .1746 
. 7576 el 
. 7269 . 0614 
.7113 . 0452 
.6791 . 2858 
. 6628 . 0643 
.6628 - „1209 
.6463 «See 
6463 - .0062 
.6295 .2332 
.6125 .3874 
.6125 .3567 
~ 5954 3778 
, 5610 . 3682 
‚5436 . 0871 
. 5262 1577 
.5262 . 1698 
. 5088 . 4586 
. 913 ‚5081 
.4913 .0971 
.4739 . 0459 
„4217 . 08/01 
. HOL6 - .0726 
. 3998 6297 
.3372 .57% 
. SERT . 3017 


. 1991 . 3129 


(u') 


. 0130 
. 0190 
. 0281 
. 0825 


. 0724 
.1223 
. 1714 
. 2660 
8851 
. 1456 
. 3590 
. 3148 
. 3033 
45551 
. 2880 
„0048 
. 1985 
. 1450 
.2391 
. 11480 
. 1226 
4424 
Amon. 
. 0266 
. 4994 
5242 
978. 
„ 2909 
dh 
. 3213 
. 7109 
. T32 
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Table Y 


. 1465 
. 1007 
.0627 
. 0545 
.0714 
.1991 
.595h 
„5782 
„5610 
. 5436 
‚5436 
‚5262 
‚5262 
‚5262 
. 5088 
. 4913 
. 4913 
. 739 
. 739 
4564 
. àoh6 
.4046 
«Seno 
30” 
. 33m 
. 3208 
„2730 
.6954 
.6125 
. 5262 
Mio, 
. 1590 
. 11465 
. 0398 


* 


9 13 1 d» Ê 161073 4 | 8 E OS TT 4 ÁÀ| 7 Voy» 5 Bo —! 2 y 0! 


cont'd) 


. 2496 
. 2023 
. 7678 
. 0825 


.2089 


.6865 
1968 
. 4165 
‚2642 
. 441448 


„261 


. 1852 


.4721 
„2441 


4452 


. 3840 


. 4361 
2191 


4511 


. 4341 
. 2886 


. 5844 


.5638 
.6438 
.5681 


63h0 
. (186 


147" 
e 3305 
. 2981 


.4759 


.7852 
. 7736 


.8839 


.8163 
.8763 
5376 
. 9420 


.4124 
.3539 
.0660 
. 3506 
1087 
. 1636 
. 4362 
. 0413 
- 4062 
. 2016 
. 3338 
. 2620 
. 0823 
s2 1O 
„3271 
. 5207 
.1136 
«2393 
. 0564 
. 3313 
.2434 
20.37 
. 2664 
.1879 
.3975 
. 3854 
.3014 
. 360T 
3152 


+ 1.1 r 2.2! ı ı 9 za+ +ı - 3 BE Fo zı ı I SEE Ố. REDE 
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Table XI - 


Least Squares Solution - 


T 


Flattening and Correction 


to Gravity Formula from Gravity Anomalies 


Forward Solution:- 


x! 
x! +67 
E 
y' 
(1) 
m1 
(2) 
u' 
(3) 
1l 
[vv] 


y 


+ 33.20 
-0.4955 


+ 19.76 
- 16.15 


+ 3.31 
- 1 


Back Solutions- 


u"' Check 


z! 


z! 


- 4,17 
+0. 0622 


- 1.29 
+ 2.07 


+ 0.78 
-0.2356 


+ 10.73 
- NOSOP 
- 0.18 


+ 10.29 
- if 


Check 


u' 


* 13.17 
-0.1966 


+ 4.87 
- 6.53 


- 1.66 
+0.5015 


+ 0.94 
+ 0.82 
+ 0.39 


+ 2.15 
-0.2089 


+ 9.65 
- 2.59 
- 0.83 
- 0.45 


+ 5.78 
- 


y 


1 


=< ess 
+ 15.41'/9 


- 460.99 


$ 51755 


50.86 
15.3656 


108853 
64.25 
- MIOS 


259.76 
25.2439 


1 


St 
203.09 
25.51 
54.26 


+ FS 


9945] 
9.4308 


+ 


+52 , 563 

-15,926.69 
- 781.49 
- 6557.36 
= SLOT 


+20, 703.39 


Check 


2 


923.80 


+ 13.7881 


HOW. 45 


+ 257.74 


+++ + 4 


+ 


22.23 
16.0997 


2m 32 
57.46 
1255 


247.33 
24.0350 


308. 74 
181.62 
26.72 
51.67 


418.73 
8.4308 


+50,548.11 
- 14,243.00 


818.83 
6243.32 
459.56 


+28, 783. 34 


x Check 


1 +9.4308 +8.4308 425.2139 424.0350 -15.3656 -16.099T7 +15.4179 +13.7881 
- 1.9701 - 1.7612 + 4.7295 + 4.2280 - 1.8541 - 1.6575 
- 5.4833 - 5.2477 + 1.4476 + 1.3854 


u! 
z! 
1 


+ /.9872 + 8.4827 


E +9.4308 48.1308 423.2738 422.2738 -16.1194 -17.1194 122.9986 +21.9987 
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69g02*0- 
GT "3 + 


GTOG *O+ 
99°T - 


9961 *0- 
LT*ET+ 
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62 *OT+ 


9G£2°O- 
9L°O + 


2200 *0+ 
LU - 


Z 


suoryenby e49T91J02 





(£) 


n 
(z) 
Z 
T a 
TE°E + (T) 
x 
cG6†t*O= Ls 
OZ°£E + L9+ x 
K X 





- IIX 3TABL 





Solution of Correlate Numbers:- 
5.78 Q,, = - 0.4825 
Qa, = 7 0-08478 
10.29 Q = + 2.12 (- .08478) = 0.1789 
aT s (0.1789 + 0.1823) # 10.29 = 0.03510 
SL da 0.78 (.03510) - 1.66 (- .08478) = - 0.4955 
Oy = (- 0.4955 - 0.0274 - 0.1407) # 3.31 = - 0.20048 
67 nh. (- .20048) - 4,17 (.03510) + 13.17 (- .08478) = 1 
Qe = (1 + 6.6559 + 0.1464 + 1.1166) + 67 = 0.13312 
5.78 Ay 7 0.3307 
qu 0.09528 
10.29 Vz + 2.15 (.09528) = - 0.2356 
Az = (- 0.2356 - 0.2019) + 10.29 = - 0.0428] 
3.31 Q y * 0.78 (- .04281) - 1.66 (.09528) = 1 
O yy = (1 + 0.0334 + 0.1582) + 3.3l = 0.36000 
5.78 Q,,, = - 0.2089 
E 0.03614 
10.29 Q,, + 2.15 (- .03614) = 1 
zz 


SITO dem = 1 
= 9711301 
Sụu 


il 


(1 + 0.0777) + 10:29 - om 73 





TT 


Standard Error of an observation of unit weight:- 


, eee 2g > = J 1456.87 = + 21.37 mgal 


n-u 
Standard Errors of the unknowns after adjustment: - 


m. = W e s cuc: 


m ' = 21.37 J 0.1331 = + 7.78 (mgal.) 

mt = 21.37 0.3600 = + 12.82 (mgal. + 1.022 x 107?) 
m ' = 21.37 0.107 = + 6.92 (mgal. & 1.022 x 1077) 
m ' = 21.37 4 0.1730 - + 8.89 (mgel. + 1.022 x 1077) 


The values for the unknowns, solved in Table XI, and their 


standard errors, are here converted to zals.:- 


æ = x’ x 107? + 0.0230 + 0.0078 gals. 
6 


6 


= y' x (1.022) x 10. 


u 


- 0.000016 + 0.000013 gals. 


<< 
| 


z' x (300889 "x 90m" + 0.0000238+ 0.0000071 gals. 


N 
i 


= ut x (1.022) x 10° =+ 0.0000096+ 0.0000091 gals. 


li 


= 
| 


The components of the longitude term are:- 


2 2 2 6 


r = 27 + W = 0.0006586 x 10 , r = 0.0000257 
u _ 9.1308 _ 
tan 2 A, = rn A 0.4052 
2X 2225, x» =12° 
O O 


Applying the computed corrections: - 


978.03C +x = 978.053 gals. 


Ya 


n 
II 


B 0.005302 + y = 0.005286 


- 0.000007 (unchanged ) 


m 
| 





78 


and adding the longitude term: - 
2 0 
0.000026 cos 9 cos 2 (A - 11”) 


we obtain the corrected gravity formula: - 
Y = 978.053 [1 + 0.005286 sin^9 * 0.000007 sin 2 Ọ 


+ 0.000026 cos^o sos 2 (A. - 399]. 


Next, the flattening f is solved from the corrected value of 6. 
At this point it should be pointed out that if the principal problem 
is to find the flattening, it would not be necessary to include 
solution of the longitude term. With slight difference the same 
result for the correction to B is obtained when only two unknowns exist 
in the normal equations. Witness, using the first reduced normal 


equation in Table XI:- 


3.31 y! + 50.86 = O 
y" = - 50.038 3. IO A 
x = -15.4 x 1.022 x 10° = - 0.000016 
B = 0.005302 - 0.000016 = 0.005286 


Now using Clairaut's Formula, (5-1), we obtain f from the new 


value of B. 


0900298955 = 2m -f- m T 

- ara "e sec! 
m e m - 86164 

m 2 T ye 637,820,000 cm. 
mo èë = \B6i6u 978.05 gal. 


39.478 x.637.02 - “j0. 2 
sior. x 918.05 ` 7,261,12À = O.0033670 





TP 


0.0052855 = 2 (0.0034678) - T [3 + zt (0. 0034618 )] 


0.008670 - f (1 + 0.0042109) 


£ _ 0.008670 - 0.005286 
- 1.00h2109 


. 0.003300 _ 
a " 0.0033698 + 0.0000083 


Wet 296.75 + 0.74 


ii 


(B = 0.0052855 + 0.0000130 


m 4 


t 7t59559 (o. 0000151) = + O.0000083) 


Of course, the result reflects the fact that only a few values 
of Ag were used, and at scattered locations. The solution by 
Heiskanen, referred to in Chapter 5, shows that the use of a great many 
points results in a value of f equal to that obtained from Helmert's 
original formula. The results in the computation given here are sur- 
prising considering the fact that the points chosen were completely 


random with no thought given to even approximately even distribution. 


7.4 Computation of Gravimetric Deflections 
7.41 Procedure and General Information. 


The procedure used for computation of the gravimetric deflect- 
ions at the two points chosen, is that involving the use of the Rice 
Circle-Rings. The extent of readily available gravity material 
dictated that the area covered in both cases was to be wiinir the 


circle-ring no. 49, having & radius of 399.0 km. 
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A template was prepared to the scale of the available anomaly 
maps, having circles of radii listed in Table I and with 10° sectors. 
The usual manner of orienting the template, which is to have zero 
azimuth bisect one of the sectors making mean azimuths of all sectors 
multiples of 10°, was not followed here because the writer's method of 
assuring proper sign determination dictated that each compartment have 
its own quadrant designation. Therefore the azimuth lines were 
numbered in multiples of 19% and the mean sector azimuths are multi- 
ples of 59, 2 

For the point in the northern hemisphere, zero azimuth is nortsu, 
and for the one in the southern hemisphere it is south. Following the 
usual practice for points in Europe, A and were considered positive 
eastward, so for both points, the azimuths were measured eastward from 
zero; thus they are measured clockwise in the northern hemisphere, and 
counterclockwise in the southern hemisphere. This will preserve signs 
automatically, using the same procedures for both. Figure 12 illu- 
strates the azimuth convention used. 

Remembering that the plumb line (below the surface) is pulled 
toward the mass surplus and thus in the direction of the more positive 
of the anomalies on opposite sides of the computation point, we can 
devise a system of assuring correct sign. Figure 13 illustrates this. 
€ is considered positive toward the nearest pole, n is considered 


positive eastward. Table XIII gives criteri& to be followed. 


.- a ë a«— 
E 


d e T alf «es JA $8 
—— A ze 


Auc M - 2 0 


-— (db — 
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HEMISPHERE HEMISPHERE 


Fig. 12 


Azimuth Convention for Computing 
Gravimetric  Deflections for the 
Points Computed in this Thesis. 
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Nearest Pole 
—> 


(AG, -AGz, ) 20 
(Ay - Agg)« o 
Quad D 
-p)>o 
Nearest Pole 
——— >> 
Ag<o 
Quad. IJE 


(A2;-A2„)<© 


(Agr - A75) >0 
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Table XIII - Criteria for Determining Effective Azimuth 


and Sign of Deflection Deflection 


[ ` 
ine sine 
Wffective s 
| Quadrants Anomaly Sradient en ` ¬ 





T and TTT Ar 
ÂrTT « O < 0 
| II and IV Ary TO ^50 
Ary >Q «TO 


Computation proceeds, beginning in the first quadrant with the 
effect of that quadrant and the one opposite, the third. After com- 
pleting the first and third, the second quadrant is begun, determining 
the effect of that quadrant and the one opposite it, the fourth. Thus 
the entire 360° in azimuth are covered. 

In order to obtain data for a graph of the effects of areas of 
different radii on the deflection of a point, the computation was made 
grouping the effects of four zones together from the point outward. 
Then accumulative effects were computed and plotted for both Ê and n. 

The inner circle, for one case 5.125 km. in radius, and for the 
other 10.15 km., was computed by the one-gradient method. Insufficient 
local gravity data at the points dictated that the three-gradient 
method would yield results of meaningless accuracy. 

The two points computed were Laiska, Finland (64° O3! N 


28° 50' E), and Wonderfontein, Union of South Africa (25° NOV S. 
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29° 53' E). The points were chosen to contribute data for the com- 
putation of the equatorial radius from the arc of the 30th meridian. 
in addition, their locations had to be such that the anomaly maps 
available to the writer would provide data to a large enough radius 
around each point to make computation by the circle-ring method 
feasible. The result was that both the points turned out to be some- 
what removed from areas of concentrated gravity stations. However, 
point-anomaly values obtained from files of the Institute of Geodesy, 
Photogrammetry, and Cartography, The Ohio State University, for the 
general area around the points, allowed the plotting of iso-anomaly 
contour maps to the scale 1" = 6 km. (1 : 240,000) for both points, and 
the writer proceeded with the computation undaunted by the lack of 
enough gravity data. 

The small scale maps used for the bulk of the computation were 


the following: - 


Laiska :  Free-air anomaly map of Finland on the scale 
1 : 1,000,000, covering the area sie 35° East longitude 
(extending into Russia) and 60° - 69° North latitude. 
Contours for the Gulf of Bothnia were sketched in from a 
smaller scale map of that area. The maps used were 
prepared from gravity determinations of the Finnish 
Geodetic Institute and were computed on the International 


Gravity Formula. Contour interval was 5 mgal. 


Wonderfontein : Bouguer anomaly lithographed map sheets of 


Union of South Africa, Northeast and Southeast sheets, on 
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scale 1 3 1,000,000, 10 mgal contours; total coverage of 
the two sheets : 22% 00' to 34° 30' S, 23° 30! to 33° 50! 
E. These maps were published in South Africa. 

Topographic maps published as World Aeronautical 
Charts by U.S. Air Force (ACIC), scale 1 : 1,000,000; 
contour interval is 1000 ft., but adequate point elevations 
to the nearest foot are given for the use to which the maps 
were put. Total coverage of sheet numbers 1275, 1276, 1299, 
1300, 1397, and 1398 Cik the same as that of the 


Bouguer anomaly maps. 


Since the South African point was estimated from a Bouguer 
anomaly map, it was necessary to use the topographic maps also, for the 
purpose of obtaining mean elevations for each compartment of the same 
Rice ring template used for the anomalies. Thus armed with the 
elevations, the Bouguer correction could be computed and subtracted 


from the Bouguer anomaly to obtain the Free Air anomaly for the com- 


partment. 
z ert 2 
^g. = Æp £g Free Air anomaly 
where ^ -  Bouguer anomaly 
8p =  Bouguer correction 


`. 7 (2 8) h 
m 


p = density of crust assumed for area 
= 2.67 
pan = mean density of earth 


= 5.576 
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Ep = = i LA (0.09406) h mgals (h in feet) 


i 


hoi 2. 03318 h mgals. 


7.42 Computation for Laiska, Finland. 


The computations for the circle-ring zones, proceeding from the 
outside zone, no. 48, and continuing inward toward the computation point 
to zone no. 27, are included in Table XIV. 

The inner circle, ro” 10.15 km., is computed below. The 
effects of both the inner De and the circle-ring zones are 


accumulated in Table XV. Curves of the effect on deflections of in-- 


creasing gravity field radius are shown in figure 15. 





Fig. 14 - Inner Circle - Laiska. 


Using the values shown in fig. 14, we obtain:- 


d£" = 0.0525 (48, - /8y) 


0.0525 (2.2 2.0) = 0.0525 (0.2) = + 8010 


dn" 


n 


0.0525 (28 - 287) 
0.0525 (- 0.3 - 4.8) = 0.0525 (- 5.1) = - 0268 
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Table XIV - Computation of Gravimetric Deflections 


at Laiska, Finland 


Ag (mgal) 
D I III I ER I T 
E 0-10 180-196 10-20 190-206 20 - M9 200-206 
48 - 1 - 30 = 5 - 30 - 20 - 20 
47 - 289 = 5 MAE - wb = IIO O 
Eg - 36 - 35 = HỘ - MP - 30 - 226 
I - III + 5 + 5 4... 
do & 07001 O. 005 O. 005 O* O10 
Br. A 185 195 25 
cos A - „9962 - „659 + .9063 
sin A - 097 - ‚2588 + .4226 
ae - 005 - 0.005 + 0:009 
an 0° 000 - 0001 + 0?00k 
h6 A - + 15 - 10 + 20 O + 20 
45 O + 25 + 5 + 20 + 5 + 20 
hl + 10 + 30 + 10 + 15 + 10 + 10 
43 + 10 + 20 + 10 + D + "HD + 10 
E + 5 + 90 + 15 + 70 + 30 + 60 
I - III - 85 - 55 - BO 
do 0.085 0.055 0. 030 
A 2 15 25 
cos + .9062 + .9659 + .9063 
sin + .0872 + .2588 + .4226 
de + .O085 + .053 + .027 
an + .007 + .014 + .013 
ho + 35 + 15 + 15 + 15 + 18 + T 
41 + 18 + 11 + 20 + 13 + 22 + 10 
LO + 18 + 10 + 20 + 10 + 25 + 14 
39 + 15 + 8 + 20 + 5 + 2h + 7 
E + 66 + yh + 75 + h3 + 89 + 38 
t- III + 9E + 32 + 51 
dO 0.022 0.032 0.051 
A 185 195 205 
cos - .9962 - .9659 - .9063 
sin - .0872 - .2588 - 26 
de - .022 - .031 - 085 
an - .000 - .008 - .082 
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sin 


an 
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a6 
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sin 


an 


48 
47 


TEL 


42 


ho 
Bp 


III 
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30 - 40 


e 
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+++ + 
~ 
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aa 
CO 
WI 


18 
22 
e 
34 


+ +++ 
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Table XIV (cont'd) 


III 
210-220 


+ 10 


+ 10 


+ 20 


10 


+++ + 
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ap 


„ 010 

215 
. 8192 
- 5736 


. 006 


+ 
O 


12 
15 
13 
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25 
. 055 
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10 
10 
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TEL 
230-240 
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O 
O 


+ SD 
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ESO, 
5736 
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++++ 
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Table XIV (cont'd) 


IV 
300-310 
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Table XIV (cont'd) 
IV II IV 
330-340 160-170 340-350 
O = 5 - 10 
- 5 + 10 - 15 
= Y + N. - 25 
+ 15 + 30 
«015 . 030 
38D 345 
+ .9063 + .9659 
- „1226 - „2588 
+ .014 + .029 
- .006 - .008 
- 5 + 15 - 10 
- 10 + 15 - 10 
- 20 + 10 - 10 
- 20 + 5 O 
- SD + 45 - 30 
+ 100 + 7 
. 100 ON 
335 345 
+ .9063 + .9659 
- 4226 - ‚2588 
+ .09 + .072 
- .Oh2 - .019 
- 20 + 7 = 10 
- 15 + 10 - 13 
> Y + 10 = 5 
- 8 + 3 - 2 
- 51 + 30 - 30 
+ 78 + 60 
.078 .060 
339 345 
+ .9063 + .9659 
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- .033 - .016 
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Table XIV (cont'd) 


II IV II IV II IV 
150-160 330-340 160-170  3ho-350 170-180 350-360 
38 + 4 - 2 + 1 O + 5 + 5 
37 + 7 - d t 3 O + 4 + 5 
36 + 10 - 3 + 4 + 2 + 1 ES 
o + ll - l * l + 2 O + 10 
5 + 32 - 10 + 9 + 4 + 10 + 28 
11 - IV + h2 + 5 = 16 
do O42 . 005 018 
A 38D 345 a7 
cos + .9063 + .9059 = .9962 
sin - .h226 - .2588 + .0872 
dt + .038 + .005 =- .018 
an - .018 - .001 + .002 
34 + 12 + 1 + 5 + 3 + 1 + 8 
33 + 15 + 1 + 10 + 3 + 3 + 6 
32 + 13 + 1 + 9 + 2 + 6 m 5 
31 + 10 O + 8 + 2 + 7 + 5 
> + 50 + 3 + 32 + 10 + 17 + 2h 
II - IV + 47 + 22 - T 
do 047 . 022 2:097 
A 339 345 175 
cos + .9063 + .%59 - .9%2 
sin - „226 ee + O02 
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30 t 9 O + 8 + 2 + 6 +15 
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Table XV 


Group Sum Accumulative Group Sum Accumulative 
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+ 0.164 
+ 0.468 
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Inner Circle 
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7.43 Computation for Wonderfontein, Union of South Africa. 


As mentioned in section 7.41, both Bouguer anomaly maps and 
topographic maps were required for this station. This situation 
prevailed from the outer zone, no. 48, through zone no. 31. The area 
within the circle described by that zone's radius, 20.09 km., was 
covered by a larger scale free-air anomaly map prepared by the writer. 

The inner circle for this station had a radius mm 5.125 km. 
Figure 16, below is a reproduction of the inner circle. Using the 


values shown, we obtain:= 


dg" = 0.0525 (Ag, - Ag,) 
= 0.0525 (99-0 - 88.5) = 020525 (10.594 0-558 
an" = 0.0525 (Ag, - 48) 


0.0525 (88.2 - 90.2) = 0.0525 (-2.0) =-0105 





Fig. 16 - Inner Circle - Wonderfontein. 
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The effects of both the inner circle and the circle-ring zones 
are accumulated in Table XVI. Curves of the effect on deflections of 


increasing gravity field radius are shown in fig. 17. 


7.44 Comments on the Results of Deflection Computation. 


The results shown in Tables XV and XVI and in Figures 15 and 
+”, are the computed effects of the gravity field around each point 
out to a radius of 399.0 km. The effect of the rest of the world, 
beyond 399.0 km. should not be ignored in practice. Time and avail- 
ability of data do not permit careful computation of this effect for 
purposes of this thesis. However, an idea of its magnitude can be 
gotten from publications containing results of investigations into 
this problem by various authors. For instance, Kaula (19), (20), 
determined the RMS uncertainties in the total deflection as computed 
from the Vening Meinesz formulas, first assuming perfect knowledge of 
gravity within given radii and none beyond, and second, ideal 
distribution of given numbers of gravity stations out to given radii. 
Taking Kaula's figures for the first case, assuming perfect knowledge 
within a radius of 200 km. and none beyond, the RMS uncertainty is 
h"1: for 500 Bm. it is 3735; for 1000 km. it is 216: ete. It should 
be noted that Kaula's computations were based on information provided 
by Hirvonen (39). 

Interpolating for 399 km. radius, an RMS uncertainty of * 376 
is obtained if Kaula's information is assumed to apply. The probable 


error is then + 0.6745 (3:6) = + 2°:h, 
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Table XVI - Summation of Effects of Inner Circle 


and Zones 23 - 48 for Wonderfontein. 


Zone o 2u Accumulative Group Sum Accumulative Ene 
Group de ae" an" an" (km. ) 

Inner + O.551 * 0.551 - 0.105 - 0.105 5.125 
23-26 + 0.515 + 1.066 + 0.013 - 0.092 1025 
27 -30 +0905 + 17972 + 0.212 « 0.120 20.09 
Bl BN + 0.128 + 2.099 + 1.099 + 1.219 39.67 
35-38 + 0.175 + 2.274 + 1.448 + 2.667 TOT 
39 - 42  - 1.88h + 0.390 + 1.139 + 3.806 151.9 
43 =- 46 - 1.919 - 1.529 + 1.222 + 5.028 291.2 
47 - 487 - 0.75% - 2.283 - 0.493 + 4.535 399.0 








5 E = = 2.283 n_= + 4.535 
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Estimates of the probable errors of the gravimetric deflections 
as computed from the data given within the field of radius 399.0 km. 
must be made, since perfect knowledge of gravity within the field can 
not be assumed. The method to be used for this estimate is that of 
the U.S.C. & G.S. and is described in Rice's paper on the 16 points 
(26). According to this method, the probable error in either component 


€ or ù, is determined from the following expression: - 


Bc = 18 E (0001 E.) 
where E, = Probable error in determining the mean anomaly 
in any compartment in a sector in mgals. 
E, = Probable error in & or n in seconds 


In making use of the above equation for this problem, the 
probable error was estimated for each zone group computed. The 
accumulated results, combined with the effect of the inner circle, each 
component of which is assumed to be no more certain than its value, and 
with the probable error for the area beyond the radius 399.0 km., give 


the following (see tables XVII and XVIII):- 


For Total Deflection O9 ;- 


Laiska E + 2.42 


t 


ii 


Wonderfontein E, + etch. 


Thus it appears that for computations based upon the gravity 
data given, and within & field of radius 399.0 km., the probable errors 
for the gravity field &re almost negligible compared with those re- 


sulting from the neglect of gravity beyond that radius. 
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Table XVII - Computation of Probable Errors in Gravimetric 
Deflections as computed for the Gravity 


Field within Radius 399.0 km. 


dy. 2 (3) (09 (5) (6) (7) . (8) 


i, 2 
Đến E, E, 0.001 E, (0.001 En) 2 (0. VOLE,,) 
Zones ri (mgal) (msal)  (r„a1) for Group 
31 (3)+(4)] 0.0015) (6) (2) (7) 
LATSKA 
27-30 h f 1 0.001 0.000001 0.000004 
31-34 4 l 7001 . 000001 . 000004 
35-38 4 l -001 » 000001 . 000004 
39-42 4 5 .005 . 000025 . 000100 
43-146 4 5 . 005 . 000085 - 000100 
47-48 2 5 . 005 . 000025 . 000050 
WONDERFONTEIN 
23-26 4 1 3 2 0.002 0.000004 0.000016 
27-30 4 l 1} 3 2 „ 0 .00000^  .OOOO16 
31-34 y 1 3 2 . 002 . 000004 . 000016 
35-38 4 ue (T 6 . 006 . 000036 - 000144 
39-42 y 5 17 11 Song . 000121 . 000484 
43-146 4 5 «3f 11 11 . 000121 . 000484 
47-48 2 5 17 11 2031 . 000121 .000242 
1 100 ft. @ 0.03378 mgal / ft. 
2 200 ft. @ 0.03378 mgal / ft. 
3 500 ft..G 0.03378«mgel / ft. 
Inner Circles:- 
Laiska Wonderfontein 
E E, - dt * 07010 + 0.551 
B 0.000100 0.303601 
n s E = an - 05208 - 07105 
E^ 0.071824 0.011025 
Bs : E e de + an? 0.071924 0. 311626 





Table XVII (cont'd) 


2 2 
be 18 s(O.001 Er)” ; E, 
(9) (10) 
Z E E 
ones t 
18 (8) 
LATSKA 
21-30 0.0000/2 
31-34 . 000072 
35-38 . 000072 
39-112 001800 
43-46 .001800 
47-48 » 000900 
For Total Deflection 6 


WONDERFONTEIN 

23-26 0.000288 
21-30 o 000288 
31-34 . 000288 
35-38 . 002592 
39-42 .008712 
43-46 » 008712 
47-48 . 004356 
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= pee. in either component 
(€ or n) of the deflection 


o. 


For Total Deflection 09 : 


(11) 
Accumulated 
L 2 


me) 
Accumulated 
E 


° v (12) 


0.000072 + 0.000 
o 000144 . 012 
. 000216 -015 
. 002016 . 0lt5 
. 003816 . 062 
. 004716 . 068 
2 
E, = 2(0. 004716 ) 
E, = 0.097 
0.000288 + 0.017 
. 000576 .024 
. 000864 .029 
- 003456 .059 
. 012168 .110 
020880 Ln 
. 025236 - 159 
E, = 2(0.025236) 
= 0.050472 
E, = 0225 
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Table XVIII - Probable Errors in Gravimetric Deflections 9 


for Combined Zones 


Zones EC a Accumulated 

Es Es 

LATSKA 

Inner Circle ức = 10.25 km.) 0.0719 0.0719 9268 

27-48 (10.25 - 399.0 km.) ` 0.009h 0.0813 0285 

Beyond 399.0 km. : (214) 5.76 5.8413 212 

WONDERFONTE IN 

Inner Circle (x, = 5.125 km.) 0.3146 0.3146 07561 

23-26 (5.125 - 399.0 km.) 0.0505 0. 3651 07604 


Beyond 399.0 km. : (2k)” 5.76 6.1251 2" 
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1:5 Approximate Correction to Finland - South Africa Arc from 


Gravimetric Deflections 


Making use of the values obtained for the gravimetric deflect- 
ions of Laiska and Wonderfontein in section 7.4, disregarding the hiệp, 
certainty arising from the fact that computation was not carried out 
beyond 399 km., a correction to the length of the meridian arc between 
the points can be made. From this, & rough idea of the corresponding 
correction to the equatorial radius can be learned. 

Geodetic and astronomic data for the two points in question 
was obtained from the Army Map Service, and was part of the data used 
in that agency's computation of "A New Determination of the Figure of 
the Earth from Arcs" (See Chovitz and Fischer (4)). Unfortunately the 
data carries a security classification for values having accuracy 
better than one minute of arc. Because of this, the actual computation 
figures cannot be published here. However, the results given reflect 
the data furnished to the writer by AMS without violating security 
regulations. 

Using latitude functions tables (either (34) or (36))for the 
International Ellipsoid, upon which the given data was based, the 
length of the meridian arc between Laiska and Wonderfontein for the 
geodetic latitudes and for the latitudes corrected by the free-air 
gravimetric and astro-geodetic deflections,were computed (See figure 
18). In effect, the two figures thus obtained were (1) the measured 
length of the arc, based on the International Ellipsoid, and (2) the 


corrected length of the arc based on astro-gravimetric data. The 





EN 


following figures give pertinent details:- 


Station Latitude Meridian Arc 
Laiska p = 6k” O03'( )" N M 
Wonderfontein @ = 25” 48r( yo Sg M, 

IND - 899 5a Me 


(M, + M_) = 9,961,448 m. 


Laiska 9, = 64° 03'( )" N (me 
= O 1 " 
Wonderfontein Đạo = E PY S (M) 
m^ = 89 EE 


(M, + M,). = 9,961,534 m. 


where Mr. and M; refer to the meridian distance from the equator to the 
point in question, q is the geodetic latitude, and Pe the gravimetric- 
ally corrected latitude. Thus "n = @+ (@t - @) ~ E 

From the above figures, we see that the total correction for 
the arc is + 86 meters. Assuming the arc to be a meridian arc of a 
sphere the size of the earth, we can thus obtain a rough correction to 


the sphere's radius:- 
ag 
SR= — x 86 
Ap 


51.3 


= A x 86 = + 55 meters = Sa 


8985 
Admittedly, this is only a rough approximation, but is close 
enough to the exact value for the data used, for illustrative purposes. 


Besides this, the accuracy of the gravimetric deflections is of doubt- 
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ful value, because of the limited gravity field used. 

Any person who is lucky enough to be in possession of data 
such as is used here, but for many points along an arc of this length, 
would proceed along the lines of section 7.2, and would doubtless 


obtain results of some significance. 


NORTH 


POLE 





LAISKA 


MERIDIAN 


WONDERFONTEIN 


Pen 19 


The Finland - South Africa Arc 





Chapter 8 - Closing Remarks 


It would be presumptuous on the part of the writer to title 
this chapter "Conclusions", because there can hardly be any conclusion 
to something that has not, in a manner of speaking, even begun. The 
lack of readily available data for pursuing & project of such 
magnitude as determining the best dimensions for the earth ellipsoid, 
prevents a significant contribution to the science, even if the ability 
to put such data to proper use were present. 

What conclusions can be made, must be based on what was al- 
ready known beforehand, and indeed could not be made without such 
knowledge, because of the nature and amount of the data used. The 
first two problems solved in Chapter 7 are what are referred to here. 
The use of only astro-geodetic deflections of the vertical is admitted- 
ly an incomplete beginning to a solution, so could hardly be expected 
to yield as good results as with the use of the gravimetric deflections. 
It is unfortunate that gravimetric deflections were not available for 
all the points used for the arc of the 98th meridian. A direct com- 
parison of methods could then have been made. It would have been 
almost as meaningless to have computed the size of the earth from only 
the astro-geodetic deflections of the 16 points in the South Central 
U.S., because of the small area covered. As it was, the computation 
by the astro-gravimetric method in the case of the 16 points, probably 
owed much of its "success" to the fact that the data given, as 


prepared by Mr. Rice, was carefully computed by him and fairly 
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accurate. 

It is hoped that the material presented in this thesis might 
provide, by its systematic presentation of the procedure to be 
followed in problems of this nature, the means whereby any interested 
reader might en the problem further, especially when in the future 
more gravimetric deflections of the vertical can be expected to be 
available. The advent of the hish speed electronic computer in 
geodetic use will hasten that day. The manual method of computing 
gravimetric deflections is so laborious that one person, in the work 
of preparation of a thesis, could not be expected to produce the 
required computations in adeguate numbers and for those stations in 
locations favorable to solution of the problem of determining the size 
of the earth ellipsoid. It is believed that the example deflection 
computations included in this thesis demonstrate this clearly. 

It is regretted that more than mere mention of two important 
methods used for solution of the problem at hand could not have been 
made in this thesis. One is the determination of the figure of the 
earth from geoid heights, and the other is the determination of the 
flattening from observations of the artificial earth satellites. 

Time, more than space, prevented their inclusion, especially since the 


writer tends to be somewhat verbose in his explanations. 
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